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ABSTRACT 

We present models that predict spectra of old- and intermediate-aged stellar pop¬ 
ulations at 2.51 A (FWHM) with varying [a/Fe] abundance. The models are based on 
the MILES library and on corrections from theoretical stellar spectra. The models 
employ recent [Mg/Fe] determinations for the MILES stars and BaSTI scaled-solar 
and a-enhanced isochrones. We compute models for a suite of IME shapes and slopes, 
covering a wide age/metallicity range. Using BASTI, we also compute ’’base models” 
matching The Galactic abundance pattern. We confirm that the a-enhanced models 
show a flux excess with respect to the scaled-solar models blue-ward ^4500 A, which 
increases with age and metallicity. We also confirm that both [MgEe] and [MgEe]' 
indices are [a/Fe]-insensitive. We show that the sensitivity of the higher order Balmer 
lines to [a/Ee] resides in their pseudo-continua, with narrower index definitions yield¬ 
ing lower sensitivity. We confirm that the a-enhanced models yield bluer (redder) 
colours in the blue (red) spectral range. To match optical colours of massive galaxies 
we require both a-enhancement and a bottom-heavy IMF. The comparison of Glob¬ 
ular Cluster line-strengths with our predictions match the [Mg/Fe] determinations 
from their individual stars. We obtain good fits to both full spectra and indices of 
galaxies with varying [a/Fej. Using thousands of SDSS galaxy spectra we obtain a 
linear relation between a proxy for the abundance, [ZMg/ZFe]ss(BaSTi)i using solely 
scaled-solar models and the [Mg/Fe] derived with models with varying abundance 
([Mg/Fe]= 0.59[ZMg/ZFe]ss(BaSTi))- Finally we provide a user-friendly, web-based fa¬ 
cility, which allows composite populations with varying IMF and [a/Fe]. 

Key words: galaxies: abundances - galaxies: elliptical and lenticular,cD - galaxies: 
stellar content - globular clusters: general 


1 INTRODUCTION 

The integrated light emitted by stellar systems such as star 
clusters and galaxies provide us with the necessary infor¬ 
mation to understand their fundamental properties, and 
how these properties arise and evolve. From observables like 
colours, line-strength indices or full spectral energy distri¬ 
butions (SEDs) we can derive the Star Formation History 
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(SFH), mean age, metal content, abundance pattern, stellar 
Initial Mass Function (IMF) and dust properties of these 
systems. 

For a quantitative study of their properties we require 
stellar population synthesis models. The method consists in 
comp aring the obse rvables to the predictions of these models 
fe.g.. iTinslevlllQ^ ). Stellar population synthesis models are 
also employed in other type of studies, e.g., to provide the 
necessary templates for kinematic or redshift measurement, 
or to determine galaxy stellar masses via mass-to-light ra- 
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tios (M/L), just to mention two applications. Evolutionary 
population synthesis models, which are regarded nowadays 
as standard, take into account the contributions of all 
possible star s, in proportio ns prescribed by stellar evolution 
theory (e.g.. IrinslevlIlQS^ l . These models combine at least 
three main ingredients that determine the quality of the 
predictions: a prescription for the IMF, a set of stellar evo¬ 
lutionary isochrones and stellar spectral libraries to predict 
a variety of observables in several bands s uch as fluxes , 
colours, M/L, or l o w res ol ution spectra (e.g., Bruzr^ 198^ 
Arimoto fc Yoshij 19861: iGuiderdoni fc Rocca-Volmerangel 


19871: 


Bressan. Chios! fc Fagotto 


19941: 


I994I: 


WorUie^ 


Ritze-Von Alyenglet^ fc Gerh a^ 

]j 994; |Vazdekis_et_alJ 199d; Fioc fc Rocca-Volmerange 
199' 


Leitherer et al.l Il999l: Maraston 20051 : IConrov fc GunnI 


Kodama fc Arimotol Il997; iMaraston 1^9^ 


2 OI 0 I) and also surface brightness fluctuations (SBFs) (e.g., 


Liu, Chariot fc Graham| |2000|: [Blakeslee, Vazdekis fc Aj hail 

200ll : lLiu. Graham fc Chariot 1^021 : Cantiello et al.ll2003l L 

In the last two decades important advances have 
been made in the field with the aid of models that pre¬ 
dict absorption line-strength indices for si ngle-age, single - 
metallicity stellar populations (SSPs) (e.g., [WorthedIl994l l 
and more complex stellar populations , built as composi¬ 
tions of SSPs fe.g.. IVazdekis et al.lll99m . These models em¬ 
ploy fitting functions that relate measured line-strength in¬ 
dices to the atmospheric parameters (Tefi, logg, [Fe/H]) 
of the library stars. The most widely used fitting func¬ 
tions are those obtained with the Lick/IDS stellar li¬ 
brary (Burstein et al.l 19841: ^ber et af] 19851: Gorgas et al.l 


brary llmjrstenL^L-aiJLLa^; |tM)er_et_aiJlia^; lijarga§_£t_aJJ 
ll99,4lWorthev et aI.lll99J : IWorthev fc Ottavianilll997l L The 

Lick/IDS system comprises definitions of 25 absorption 
line-strengths in the optical spectral range with varying. 


ting functions in the optical range (e.g., 

Buzzoni 

19951: 

Gorgas et al.lll999liTantalo. Ghiosi fc Piovan 

120071 Sr.hiavonI 

200711 and in other spectral ranges fe.g..lGenarro et al.l 

2 OO 2 I; 

Marmol-Oueralto et al.ll2008l: iGenarro et al.l2009h. The ma- 


jority of these model predictions are implicitly based on a 
scaled-solaiE] heavy element mixture around solar metallic- 
ity. This is due to the fact that these libraries are composed 
of stars that are imprinted with the abundance pattern of 
the solar neighborhood, as a result of their experienced SFH. 
Such types of predic tions have also been developed in other 
spectral ranges (e.g. , ISaglia et al.ll2002l : IVazdekis et al.l[2004 : 
ICenarro et al.ll2f)09l L 

There are a variety of approaches that have been em¬ 
ployed in the literature to extract the information con¬ 
tained in spectral line indices. The most popular method 
is to build key diagnostic model grids by plotting an age- 


fe.g.. Mgfo, (FelL to estimate these parameters fe.g.. iR.ose 

1994 

IWorthevlll994 

Vazdekis et al.lll996l:lKuntschner et al. 

2 ooe 

SchiavonI 2007 

), including estimations of the abun- 


iThomas. Maraston fc Bendeij l2003l : ISchiavon 2007 1. There 


are other methods that employ as many Lick indices as 
possible and simultaneously fit them in a sense (e.g., 


^ Some authors use the term “solar-scaled” as opposed to scaled- 
solar, and we regard these terms as interchangeable. 


IVazdekis et al.l Il997l : IProctor. Forbes fc Beasl^ l2004l l. An 
alternative method is the use of Principa l Component Anal¬ 
ysis (e.g., ICovino. Galletti fc Pasine^m 1 19951 : IWild et ^ 

I 2 OO 9 II . 

Unlike colours, these in dices are insensit ive to the 
effects of dust extinction fe.g.. l^acArthurll2005l l and allow 
us to study the age (as a proxy for the SFH) and metallicity, 
[M/H] (i.e., the total metal content). In addition they also 
allow us to study the abundances of various ele ments 


rrom ine strenfftns oi tnese aosor 
19851. I 1994 Sanchez-Blazauez et al. 

)tion lines le.e.. nose 

2 OO 3 I: Garretero et al. 

2004 iTrager et al. 


19981. l200nl: IProctor fc SansomI 120021: 

Thomas et alj 2005 

SchiavonI 20071: Gonrov & van Dokkum 


2012bh . In fact. one of the most important results has 


been the finding that massive early-type g alaxies (ETGs 
show an overabundance of JMg/Fe] (e.g., IPeletierl Il989l : 
IWorthev. Faber fc Gonzaled ll992l L Such departure from 
the solar abundance pattern is generally attributed to short 
formation time-scales (<lGyr) for the bulk of their stellar 
populat i ons (see the reviews of e.g., iTrager et al1ll998l.l200nl: 
iRenzinil 12004 : IConrovI l2018l . and references therein). The 
Magnesium, an a-element, is ejected into the ISM by type-II 
supernovae resulting from massive stars on time-scales of 
~10Myr, whereas most of the iron is released on time-scales 
of Gyr, when SNIa explosions of intermediate-mass stars 
in binary systems occur. Although other a-elements are ex¬ 
pected to follow magnesium there are except ions to the rule 


sucn as caicium, wnic 
1997l:IWorthevll 19981: 

n nearly tracKS re le.e.. vazaeKis et ai. 
Saglia et al.ll2002l: IGenarro et al. 20031: 

Falcon-Barroso et al. 

l2003al. Thomas. Maraston & Bender 

2003bl: Yamada et al.l 2006: Worthev. Ingermann fc Serven 

2 OIII: iGonrov fc van Dok 

■cum 2012bl: 

Johansson. Thomas & MarastonI 

2012llLa Barbera et al. 


20131^ These abundances and the abundance ratios 
of other elements can therefore be used to inves- 
tigate the starformation timesca les of ETGs (e.g., 
IConrov. Graves fc van Dokkum|[2^L 

Aiming at obtaining more quantitative determinations, 
models with varying abundance ratios were developed 
with the aid of theoretical stellar spec tra (|TriDicco fc Belli 
I 1995 I : iKorn. Maraston fc ThomasI 1200,4 1. The responses ob¬ 
tained for the Lick indices to variations of the indi¬ 
vidual element abundances, i.e. the so called response 
functions, were used to correct the mo del predictions 
based on the empirical stellar lib r aries ( e.g., Vaadekiset^^jJ 


W9ijTantalOj_^hiosi& 


LlTVager et^ 


Thutalo_&_Chiosi [20041 : iThomas et al.lf2003l : ISchiavonI 

ioi4l 


200 c 


2007 


Ghung et al. 201311. These models allowed us to achieve a 


better understanding of what these line-strengths are made 
of, as none of these rather strong atomic and molecular lines 
are contributed to by a single element. Although with the 
use of these indices it has been possible to partially lift the 
age/metallicity degeneracy affecting the colours, these in¬ 
dices are not comple tely free from this degeneracy llWorthevI 
11994 lArim(^ll99filL 

This degeneracy is due to the fact that an age increase 
causes an effect on the theoretical isochrone and on line- 
strength indices that can be mimicked - at constant age 
- by a metallicity increase. The effects of this degeneracy 
are stronger when low resolution indices are used, as the 
metallicity lines appear blended. For such a reason alter¬ 
native indices - that were thought to work at higher spec- 
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tral resolutions, such as those o f e.g., IRos i d 19851. I1994| : 
IServen. Worthev &: Brile-y^ (120051 1 have been designed to 
alleviate the problem. 

New means for improving the stellar population anal¬ 
ysis have been achieved by predicting full spectral energy 
distributions (SEDs) at moderately high resolution, rather 
than just a number of indices. Such modelling has been 
possible with the development of extensive empirical stellar 
spectral libraries with flux-calibrated spectral response and 
good atmospheric parameters coverage. Among the most 
widely employed st ellar librari es in the optical range ar e 


widely employed st ellar librari es m tne optical range ar e 
those of l.Tonea (Il999l'l. ELODIE (Prugniel fc Soub iranl2nnill. 
STEL IB l|Le Borgne et al.l l200,8l'l. Indo-US dValdes et al.l 
|2004 1 and MILES dSaiichez-Bjazq uez et aid I 2 OO 6 I I and 
more recently XSL dChen et al.ll201 4h. Mo d els that employ 
such li braries are, e.g.. IVazdekisI lll999h: ISchiavon et al.l 


d2002tl: Bruzual fc CharlotJ d2003li: LeJ3orgng_eL^ 


^ iMaraston fc StrombackI 1 2011 
Conroy fc van DokkumI d2012al 'l . There are also 
pirical librari es i n othe r spec tral ranges such as, e.g. 


ICenarro et al.| d2001al 3); ICiishing. Ravnier fc Vaccal 


l|2005h : Ravnier. Cushing fc Vaccal d200 9l'). that are also 


implemented in thes e mod els dVazdekis et al.l l2003l : 
I Conroy fc van DokkumI l2012al l. Alternatively, theoretical 
stellar libraries at high spectral re solutions have also 
been developed for this purpose (e.g., Bj3£bia_eLal| 
Murphy fc Meiksinl l2004llzwitter. Castelli fc Munaril 


Coelho et al.l l2005l: iMartins et a l2005l : 

2 OO 5 I : iRodrfguez-Merino et al. 200^ 


Martins fc Coelh3 l2007l: 


200 . 


200 


jMunarTe^^i 


PalaciosetalJ 20101 : Kirbvl 20 111 : Ide Lavernv et all 


^__,__ jFreni^x et al 

jLeithereretay 


2 OIOI 


2 OI 2 I 


Coelhc I 2 OI 4 I'). Model s that m ake use of these libraries are. 


e.g.. i Le itherer et al.| (|l999|L ISchiavon. Barbuv fc Br uzuall 

dho et al.l d2007ll: 


d2(30Cl Gonzalez-Delgado et al. d2005ll:ICoelho 


iBuzzoni et al. d2009ll: ' bee et al.l ( 2009l l : Percival et all 

d2009ll : IMaraston fc StrombackI ( 2011 ). 

These models provide ways to alleviate the main de¬ 
generacies hampering stellar population studies. For ex¬ 
ample, with these SEDs it has been possible to define 
new indices with enhanced abilities to disentangle the 


age from the metallicity (e.g., Vazdekis fc Arimoro 19991 : 


iBruzual fc Charl^ l20o4 Cervantes fc VazdekisI 1200 91'). or 


even t he IMF (e.g., |Schiavon__eL_aL I l200d: ICenarro et al.l 


2001^ Vazdekis_e£_^ 1200^ Conro^_&raBCDokkun^ 201^ 


Vazdekis et al Jl20ld;E Barbera et alJl201.'d : ISpiniello et all 


20141) . It is becoming common practice to use these 


model SEDs with newly developed algorithms that aim 
at constraining and recovering the SFH by means of full 
spect r um-fitting approach (e . g..[Panter. H eavens fc Jimen^ 
2 OO 3 I. Cid Fern andes et ^ 2005l:~ Qcvirk et all 2006a h 


20061: Yamada et al. 

I 2 OOGI: Maraston & StrombackI I 2 OIII: 

Vazdekis et all 

2012 

ConrovI 2013) and also nhotomet- 

ric data (e.g., 

Maraston et alJ l2009l: iPeacock et alJ l201ll: 


on theoretical libraries, these models are generally re¬ 
stricted in their SSP parameters coverage, as most of the 
stars come from the solar neighbourhood, with a charac- 
teristic, and limited, elemental abundance ratio coverage 
llMilon e. Sansom fc San chez-Blazauej[20lll : see also figure 1 
in IWalcher et al. 20091 which illustrates the effect of the 
solar-neighbourhood pattern on models based on empirical 
libraries). The quality of the resulting model SEDs relies to 
a great extent on the atmospheric parameters coverage of 
the library. 

As most of the stars of these empirical libraries belong 
to the solar neighborhood, the resulting model SEDs have a 
characteristic abundance pattern: at low metallicity the pre¬ 
dicted spectra are enhanced in a-elements, whereas around 
the solar value the models are nearly scaled-solar (e.g., 
lEdvardsson et aI]ll993l : lMcWilliamlll99^ 1. In addition, most 
of these models employ scaled-solar isochrones and therefore 
they are not self-consistent at low metallicities. These mod¬ 
els, which can also be regarded as ’’base models”, present im¬ 
portant limitations for studying massive ETGs (or any other 
system that deviates strongly from the chemical pattern of 
the solar neighbourhood), which reside in a high metallicity 
regime and show enhanced [Mg/Fe] abundance. It has been 
show n that by using a-enhanced stellar evolution models 
(e.g., ISalaris. Groenewegen fc Weissl|200C ) and the stellar 
theoretical spectra (e.g., iTripicco fc Belli 19951: iKorn et al.l 


I 2 OO 5 I : ICoelho et al.ll2007tl . line indices such as Mgfe and (Fe) 
are affected. However, it is worth noting that the base mod¬ 
els have been successfully applied to these galaxies. In fact 
it is possible to obtain a good proxy for the [Mg/Fe] abun- 


vsis (IVazdekis et al. 

l2001al: ISanchez-Blazauez et al.l l2006bl: 

Yamada et al.l 

2006; 

de la Rosa et al.l 

2007: Michielsen et al.l 

20081: Vazdekis et al. 2010l: La Barbera et all 20131). This 


proxy leads to a linear relation with the abundance ra¬ 
tio estimated with the aid of models that specifically 
take into account such non-solar elemen t mixtures (e.g., 
T.Mhalo_et_aljJjj9j;_ [Tl^^ all l2003l : iLee fc WorthevI 


2 OO 5 I: ~aves fc Schiavonll2008l). 


There are some examples of models predicting SSP 
spectra at moderately high resolution with varying abun¬ 
dance ratios. Spectral SSP models that include varying 
abundance ratios both at the isochrones a nd st ell ar spectra 


have been presented in ICoeUio et al.l (l200^}]_ Lee^jt^ 


(2009 

): iPercival et al. 

(2007 

): Prugniel et al. 


Ch^ngarian et al. I I 2007 I : iKoleva et al.l l2009l : iToieiro et al’l IConrov fc van DokkumI 


2oT 3). This method represents an alternative to the more 


S Also iGervantes et al 
; IWalcher et al ' 

) present models 


1 20091) : 
with 


classical approach o f fitting a selecte d number of line- 
strength indices (e.g.. lTrager et ahll 19981 ). Furthermore these 
model SEDs have been used in a variety of applications. 
For example, these SSP SEDs have been shown to improve 
the analysis of gal axy kinematics for both absorption and 
emiss ion lines (e.g.. lFalc6n-Barroso et ^l2003bl : ISarzi et al.l 
l200d) . Models based on empirical libraries are free from 
the uncertainties in the underlying model atmospheric cal¬ 
culations and tend to provide go od fits to both, absorp- 
tion line-strengths and spectra (e.g. lSanchez-Blazauez et al.l 


varying abundances but considering only the effect of 
the stellar spectra. There are also examples where only 
a-enhanced isochrones are employed, ignoring the ef- 
fects of the stellar spectr a in t he line-strength s (e.g ., 

I Weiss. Peletier fc Matteuccil Il995l : IVazdekis et all 12003 ). 
Self-consistent models can only be achieved when both 
the spectral and the evolutionary effects are correctly 
taken into account, an d thus require the use of both 
theor etical spectra (e.g., iGoelho et al.l [ 2 OO 5 I : iMunari et al.l 
2005|), and isochrones with varying [a/Fe1 abundance (e.g., 
Salaris fc WeissI Il998l : ISalasnichl ]2()0(]l : IVandenBerg et aP 
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200^j_20H- Kim et al. 200 IPietrinferni et alJ l2006l : 
Coelho et alJ I 2 OO 7 I : Potter et alJ 20081 '). Tb our knowl¬ 


edge, the three self-consistent models for a fair range 
of ages and metal lic ities avail a ble in the literature are 
ICoelho et al.l |20o3); iLee et ahl ll2009h and IPercival et al.l 
1 200lt l. However, both these models carry a caveat: the stel¬ 
lar spectral libraries employed are theoretical, and as such 
they carry wavelength-dependent syst ematic deviations 
when compared to observed spectra llMartins &: Coelhd 
I 2 OO 7 I : |Bertone_e£^ 20M Coelhd l2014h. Inspired by th e 
work of Cervantes et al.1 "1 200711 md Prugniel et al.1 (lioOif) . 
and in order to bring together the predictive power of 
fully-the oretical models and the accuracy of the base 
models, IWalcher et al.l (l2009l l proposed a method to 
use the predictions of fully-theoretical SSP models only 
differentially. 

Here we further develop our models to predict SSP 
spectr a with varying [g/Fe] ratios based on the MILES l i- 
brary (ISanchez-Blazauez et al.ll2006l : ICenarro et al.ll2007^ 1. 
For this purpose we apply a differential correction to our 
empirical SSP models w hich is obtained from the t heoretical 
stellar spectral library of ICoelho et al.1 ||200 !tI . l2007l 'l . We also 
employ the scaled-solar and a- enhanced isochrones of BaSTI 
l|Pietrinferni et all l2004 l200d ') . This work is the second of 
the series that we started with I Vazdekis et al.l l|201(t l. here¬ 
after Paper I, where we presented our base model predictions 
that were fully based on the MILES library. As the param¬ 
eter coverage of MILES constitutes a significant improve¬ 
ment over previous empirical stellar libraries, these models 
are safely extended to intermediate-aged stellar populations 
and to lower and higher metallicities. 

The layout of the paper aims to permit the reader to fol¬ 
low the model construction from the very first steps, through 
the characterisation of the model outputs and finally to 
provide some illustrative applications. In Section we de¬ 
scribe the main model ingredients, which include the IMF, 
the isochrones and the empirical and theoretical libraries 
employed here. In Section [3] we describe our approach for 
constructing these SSP SEDs. The quality and behaviour of 
these SSP spectra are assessed in SectionlH In Section[5]and 
Section [6] we show the behaviour of the line-strength indices 
and broad-band colours as measured on the resulting SSP 
spectra, respectively. In Section [7] we compare our models to 
predictions from other authors. In Section [8] we apply these 
models to data of Galactic Globular Clusters and ETGs to 
illustrate the predictive power of these new models. We com¬ 
pare the colours, the full spectra and the line-strengths of 
ETGs with our models. We also calibrate a proxy to estimate 
the [Mg/Fe] abundance ratio, relying solely on scaled-solar 
predictions. In Section[9]we present the new features that we 
have introduced in our web-page to make it publicly avail¬ 
able these predictions, including more complex models with 
arbitrary star formation histories. Finally in Section [TOl we 
provide a summary. 


2 MODEL INGREDIENTS 

In this Section we describe the main ingredients in our new 
stellar population synthesis models, namely the IMF, the 
isochrones and the stellar spectral libraries. We summarise 
the chief aspects of these ingredients that are particularly 


relevant to this work. We also describe the details and the 
treatment that have been required for implementing these 
ingredients in the models, which include new isochrones that 
were specifically computed for this work. 


2.1 IMF 


Several I MFs are considered: t he two power-law IMFs de¬ 
scribed in I Vazdekis et all dlDhfil L i.e unimodal and bimodal, 
both characterised by the logarithmic slope, F and rJ!|, re¬ 
spectively, a s a free param eter, and the multi-part power- 
law IMF s of Krouoal l|200 j l. i.e. universal and revised. The 
ISaloeterl |l95^ IMF is obtained by adopting the unimodal 
case with slope F = 1. Although not identical, the 
Kroupa Universal IMF is very similar to a bimodal IMF 
with slope Fb = 1 .3. An extensive descr iption of these defi¬ 
nitions is given in I Vazdekis et all (l2003lf (Appendix A). We 
set the lower and upper mass-cutoff of the IMF to 0.1 and 
100 Mq, respectively. It is worth noting that the very low- 
mass (< 0.6 Mq) tapered ’’bimodal” IMF has been found 
to be consistent with the line-strengt hs and mass-to-light 
(M/L) ratios of massive ETGs (e.g., Ferreras et al.l I 2 OI 3 I : 

' 2014lL On the con¬ 


iLa Barbera et ah I I 2 OI 3 I : ISpiniello~ et al. 


trary, the single power-law ’’unimodal” IMF leads to ex¬ 
tremely high M/L values that are not supported by inde- 
pendent gravitational leasing and dyn amical estimates (e.g., 
iTVeu et'^l2000l : ISDiniello et al.ll2012ll . 


2.2 Isochrones 


We employ two sets o f theoretical isoch rones. The 
scaled-solar is ochrones of Girardi e t al.l ll2000lf (hereafter 
PadovaOO) andIPie trinferni et al.l (l 2004h and the a-enhanced 
models of IPietrinferni et al.1 l ~06ll . The latter two will be 
regarded hereafter as BaSTI. 

The Padova isochrones cover a wide range of ages, from 
0.063 to 17.8 Gyr, and six metallicity bins {Z = 0.0004, 
0.001, 0.004, 0.008, 0.019 and 0.03), where 0.019 represents 
the solar value. The isochrones also include the later stages 
of stellar evolution, using a simple synthetic prescription 
for incorporating the thermally pulsing AGB regime to the 
point of complete envelope ejection. The range of initial stel¬ 
lar masses extends from 0.15 to 7Mq. The input physics 
of the se models was updated with respect to lBertelli et al.1 
l|l994h with an improved version of the equ ation of state, 
the opacities of I Alexander fc FergusonI l|l994l 'l and a milder 
convective overshoot scheme. A helium fraction was adopted 
according to the relation: Y « 0.23 -b 2.25Z. 

We also use here the BaSTI theoretical isochronefl. The 


^ Note that we have changed our notation for the IMF slope from 
that we us ed in o ur previous papers, i.e . fi (e.g.. IVazdekis et al.l 
Il996l . l2003l .PaDer l lVazdekis et al.1 ll2012h '). We now adopt F and 
Ff, for the slopes of the unimodal and bimodal IMF shapes, re¬ 
spectively, to be consistent with recent works on the variation 
of the IMF. Note also that in the new notation we distinguish 
the slopes of these two IMF shapes as they apply to different 
functional forms. 

® In reality it would be more accurate to use F = 1.35 to match 
the Salpeter slope 

The BaSTI database is publicly available at 
http: //www.oa-teramo.inaf.it/BASTI 
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original sets of models have been supplemented by addi¬ 
tional computations specifically performed for this project 
(see below for more details). A description of the complete 
BaSTI database as well as of the adopted p hysical inputs and 
numerical a s sump t ions can b e found in f^ trinfern i et al.l 
(|2004 I 2 OO 6 I. l200i l2013|j and ICordier et al.l (1200711 . We re¬ 


fer the interested reader to these references. Here it is 
sufficient to note that the BaSTI archive is based on an 
updated physical framework and that it has been exten¬ 
sively tested with many observational constraints based on 
both resolved (mainly eclipsing binaries and cluster color- 
ma gnitude diagrams) and un resolved stellar pop ulations 
(see lPietrinferni et alJ[2003 and IPercival et ^l20od l. In this 
work we have adopted the non-canonical BaSTI models 
with the mass loss efficiency set to 77 = 0. 43 for both the 
scaled-solar and a-enhanced heavy elements mixture. The 
initial He mass fraction ranges from 0.245, for the more 
metal-poor composition, up to 0.303 for the more metal- 
rich composition, with AYjAZ « 1.4. All the metallicity 
grid points available in the archive have been used. How¬ 
ever in order to better sample the whole metallicity range 
from very metal-poor to the most metal-rich stellar popula¬ 
tions we have added an additional grid point correspond¬ 
ing to Z = 0.024, Y = 0.279 (i.e. [M/H] = -f0.15) for 
both scaled-solar and a-enhanced mixture. This metallic¬ 
ity is particu larly relevant for study ing massive early-type 
galaxies (e.g. ET Barbera et al.lEolij '). Therefore the BaSTI 
metallicity grid adopted in present analysis was based on 12 
grid points: Z = 0.0001, 0.0003, 0.0006, 0.001, 0.002, 0.004, 
0.008, 0.0100, 0.0198, 0.0240, 0.0300 and 0.0400. For each 
metallicity, the isochrone age range covers the interval from 
0.03 to 14 Gyr, with a fine age grid. 

We note that the atomic diffusion of both helium and 
metals was properly accounted for when computing the solar 
model in order to be able to match accurately the helioseis- 
mological constraints. The best match to the depth of the 
convective envelope (0.716 Rq) of the present solar envelope 
He abundance (Y= 0.244), and of the actual (Z/X) ratio 
[(Z/X= 0.0244], lead to an initial He abundance and metal¬ 
licity (Yq = 0.2734, Zq = 0.0198), which was adopted for 
the solar metallicity model set. 

We have used here the ‘AGB-extended’ version of the 
BaSTI isochrones that account for t he Asymptoti c Gian t 
Branch (AGB) stage as provided bv ICordier et al] ll2007li . 
This phase includes the entire Thermally Pulsing AGB (TP- 
AGB) regime th rough a synthetic-AG B treatment following 
the technique of llben fc TruranI ill978l l. In this approach the 
mass- loss processes are accounte d for by using the formu¬ 
lae of I Vassiliadis fc WoodI (Il993ll . The evolution with time 
of the CO core mass and surface luminosity is accounted 
for by adopting the analyti c al rel ationships provided by 
IWagenhuber fc GroenewegerJ (Il998l l. while the trend with 
time of the effectiv e temperature is pro vided according to 
the relationship of IWagenhubeil (Il996l l. We wish to em¬ 
phasise that our synthetic-AGB treatment does accounlQ 


® We note that r} rep resents the fre e parameter present in the 
Reimers mass loss law llReimer adopted in the BaSTI stel¬ 

lar model computations. 

^ We would like to emphasise that the claim provided by 
iMarigo et al.l ll2008h that the AGB-extended version of the BaSTI 
isochrones does not explicitly account for AGB nucleosynthesis, 


for the occurrence - and, hence for the related evolution¬ 
ary effects - of the third dredge-up and hot bottom burn¬ 
ing when appropriate. We notice that this simplified model 
for the TP-AGB phase has been shown to be adequate for 
matchin g several integrated p roperties in the near-IR as dis¬ 
cussed in Cordier et al.l ll2007l l. In addition, in order to assess 
the reliability of the theoretical framework adopted in the 
present work, we would like to call the attention of th e reade r 
to the recent analysis performed by ISalaris et al.l ll2014ll . 
They presented stellar population synt hesis models based 
on fu lly evolutionary AGB computations dWeiss fc Fergus^ 
I 2 OO 9 II (i.e. without the need of adopting analytical relations 
as in the synthetic-AGB treatment) and compared their 
model predictions with integrated colours from other mod- 
els. They found sig nifica nt discrep a ncies with results from 
iMarigo et al.l (|200^ and iMarastonI ll2005lf1 . whereas an ex¬ 
cellent agreement was found with the BaSTI AGB-extended 
models. 

The lower mass limit in the standard BaSTI isochrones 
is 0.5 M 0 . However, for the aim of the present study, we have 
extended the isochrones to the very low-mass (VLM) regime 
by including additional stellar models down to a lower mass 
limit of 0.1 MfT). The adopted V LM models are the same 
presented in iGassisi et al.[ ll200d '). which are based on an 
updated physical framework for the evolutionary computa¬ 
tions of VLM stellar models. We have made a significant 
effort to obtain a smooth transition when passing from the 
low mass regime (BaSTI models) to the VLM ones. This has 
been obtained by computing some additional stellar models 
in the transition mass regime (around « 0.5 — 0.6 Mq) by 
using the same physical inputs adopted for the BaSTI model 
computations but the same outer boundary co nditions used 
in th e VLM stellar model computations (see ICassisi et al.l 
I 2 OOOI for more details on this issue). Our temperatures for 
these stars are cooler than those of PadovaOO, and warmer 


i.e. the effects of the third dredge-up and hot-bottom burning, 
in the underlying AGB models, is not correct. The analytical 
relationships for the trend of the surface luminosity and effec- 
tiv e temperature as a function of the He core mass provided 
bv IWagenhuber & GroenewegerJ lll998ll and IWagenhuben lll996ll 
are adopted in the BaSTI synthetic-AGB treatment, do account 
for the occurrence of these physical processes during the AGB 
stage. In order to provide a correct descript ion of the situation , 
we wish to note that the computations by IWagenhuberl lll996ll 
were performed without accounting for low-temperature radiative 
opacities with various G/O ratios - because these opacity tables 
were not available at that time. Notwithstanding, the effect of the 
change in the envelope G/O ratio induced by the third dredge- 
up was mimicked by adopting a radiative opacity at a constant 
heavy e lements distribution, b ut al lowing the global metallic ity to 
change. [Cristallo et al.l ll2007h and IWeiss &: FergusonI ll2009l') have 
shown that, although the use of opacity tables that properly in¬ 
clude the enhancements of G is more appropriate, the predictions 
of the behaviour with time of Tgff as provided by fully consistent 
AGB models are not significantly different from those predicted 
by models that account for the opacity effects of the third dredge- 
up by adopting a variable global metallicity but constant G/O 
ratio. 

^ A comparison among these population sy nthesis models an d 
observations has been already performed bv iNoel et al.l ||2013|) . 
and it was found that these models disagree - being in general 
too red — with the integrated colours of superclusters for ages 
between 10^ and ~ 10^ years. 
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than those of Polg_et ^ (Il995l l , wh ich were implement ed in 
I Vazdekis et alj J l996h . and those of iDotter et alJ (l2008l 'l (see 
also I An et al.ll2009h . As their temperatures are higher the 
PadovaOO VLM stars are brighter than the ones we employ 
here for the BASTI models (see figure 5 in IVazdekis et alJ 

l2ni2h . 

In Fig. n we plot a set of BaSTI isochrones for two 
selected metallicities and ages. We see that the higher 
the metallicity, the greater the difference between the 
scaled-solar and the Q-enhanced models, with the scaled- 
solar isochrones being slightly cooler th an the a-enhanced 
isochrones Isee IPietrinferni et al.ll200(ih . These differences 
arise since we are comparing the two sets of isochrones at 
constant global metallicity; when comparing the two sets at 
constant iron content, the a-enhanced isochrones would ap¬ 
pear cooler that the corresponding scaled-solar isochrones as 
a consequence of their larger global metallicity and, hence, 
hi gher opacity at low temperatures, (see, e.g., figures 1 and 2 
. Coelho et al.ll20(y^ and figure 1 in lLee. Worthev fc Dott^ 


l2009bll . A mOTe elaborated discussion on these effects, ir 
eluding those when varying the abundance ratio of individ¬ 
ual g -elements can be found in iLee. Worthev fc Blakesl^ 
l|201(t l. While there are clearly differences between the two 
sets of isochrones, it is worth noting that these differences 
are always small. 

We use the theoretical parameters of these two 
sets of isochrones (Tefj, logg, [M/H]) to obtain stel¬ 
lar fluxes on the basis of empirical relations between 
colours and stellar parameters (temperature, gravity and 
metallicity), instead of using theoretical stellar fluxes. We 
m ainly use the metallicity-dependen t empirical relations 
of I Alonso. Arribas fc Martfnez-Rogerl l|l99fil . fl999l 'l: respec¬ 
tively, for dwarfs and giants. Each of these libraries are 
composed of ~ 500 stars and the temperature scales are 
based on the IR-Flux method, i.e. only marginally depen¬ 
dent on model atmospheres. We use the empiri cal compi¬ 
lation of iLeieune. Cuisinier k. Bused ll 19971 . fl9^ (and ref¬ 
erences therein) for the coolest dwarfs (TefF<4000K) and 
giants (Teff;^3500 K) for solar metallicity, and also for stars 
with temperatures above ~8000K. We use a semi-empirical 
approach for these low temperature stars for other metal¬ 
licities. For this purpose we combine these relatio n s and 
the m odel atmosphere p redictions o f Bessell et al.l (1 19891 . 


[i 993) and the library of iFluks et al.l (1994). We also em- 


ploy the metal-dependent bolometric co rrecti ons given by 
lAlonso, Arribas fc Marti'nez-Rog^ (1 19951 1 and I Alonso et all 
( I 999 II for dwarfs and giants, respectively. Finally, we a dopt 
BCq = —0.12. By assuming V© = 26.75 llHaveslll985fl . the 
solar absolute magnitude is Mvq = 4.82, therefore, for the 
adopted normalization, M^oIq is given by My© + BCvq = 
4.70. 

In the present work the integrated colors and spectral 
indices based on the BaSTI stellar models are compared to 
those obtained on the basis of the Padova models. Therefore 
we consider it worthwhile to briefly discuss some of the prin¬ 
cipal differences between the two theoretical evolutionary 
frameworks. Fig. [2] shows a comparison for three represen¬ 
tative metallicities (Z=0.001, Z=0.004 and solar) and three 
ages (10, 1 and 0.1 Gyr). A detailed comparison between the 
scaled-solar BaSTI models and the PadovaOO isochrones has 
been performed in IPietrinferni et ahl ll2004ll , and we refer 
the reader to section 6 of this paper for details. Here it is 


sufficient to note that there is a systematic offset between 
the BaSTI and Padova models in the RGB Tefi scale, with 
BaSTI being hotter models at solar metallicity, whereas at 
lower metallicities the trend is reversed and BaSTI is cooler 
than PadovaOO. There are also significant differences in the 
brightness of the isochrones (at a given age) for ages lower 
than about 2 Gyr. These differences can be fully understood 
when accounting for the differences in the adopted physical 
inputs and in the different treatment of core convective over- 
s hooting during the core H -burning stage (see the discussion 
in IPietrinferni et al.ll200j 'l. We also see important differences 
in the AGB stellar evolutionary phase, as illustrated by the 
middle panels of Fig. O These late evolutionary stages are 
brighter in the Padova models at 1 Gyr, whereas the BaSTI 
AGB extends towards cooler temperatures, with such trends 
depending on the metallicity. This difference is quite well 
known as these stars reach the tem perature regimes wher e 
the B — V colour saturates (see, e.g.. lWorthev &: L^l201ll l. 
Note that this difference is further emphasised at younger 
ages. Similarly, the two sets of models show important dif¬ 
ferences during the core He-burning stage: for old ages (i.e. 
larger than ~ 3 — 4 Gyr) this is due to both the differ¬ 
ences in the adopted physical framework and (mainly) in 
the adopted efficiency of mass loss along the RGB stage; 
while for younger ages the origin of the differences in the 
‘blue loops’ - characteristic of the core He-burning stage of 
intermediate-mass stars - could be related to many different 
contributions (physical inputs, treatment of convection, the 
details of the heavy ele ments mixtures, etc.) as discussed in 
detail bv lGassis 1 (l2004lb 

Another interesting feature that is relevant to IMF stud¬ 
ies is the difference in effective temperature between these 
two sets of models for the VLM stars, with the Padova 
models being hotter in all the panels of Fig. O This is re¬ 
lated to the different physical inputs to the models - mainly 
in the equation of state and outer boundary conditions - 
adopte d for computing th e two sets of VLM stellar mod¬ 
els (see iGassisi et al.l I2OOOI and references therein for more 
details on this topic). 


2.3 Stellar spectral libraries 

2.3.1 MILES 

The new models presented in this work, for both the scaled- 
solar and a-enhanced mixture, are computed - in the optical 
range - on the basis of the M ILES stellar spectral librarjjf] 
(ISanchez-Blazauez et al.l I2OO6II . This library was prepared 
carefully to be implemented in our models (Paper I). For a 
detailed de scription we refer the r eader to these two papers 
and also to lGenarro et al] (l2007al L for the determination of 
the stellar parameters. Here we briefly summarise the most 
relevant aspects of the library. 

The MILES library is composed of 985 stars observed at 
the INT 2.5m telescope on El Roque de Los Muchachos ob¬ 
servatory at La Palma, Spain. Two gratings were employed 
to cover the full optical range at 2.51 A (FWHM) resolution, 
which is constant with wavelength fsee lEalcon-Barroso et all 
I2OIII for a full description). The spectra have very high 


The MILES library is publicly available at http://miles.iac.es 
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Figure 1. Comparison of a-enhanced ([a/Fe] = +0.4)and scaled-solar BaSTI isochrones. Two representative metallicities Z= 0.024 (left) 
and Z= 0.001 (right) and two ages 12 Gyr (top) and 2 Gyr (bottom) are shown. 


signal-to-noise (S/N) which is usually well above 100, with 
the exception for the stars that belong to globular clusters 
such as M 71, which were included to improve the parameter 
space coverage. The spectra were carefully flux-calibrated; 
each individual star was also observed through a wide slit 
in order to avoid selective flux losses due to differential re¬ 
fraction effects, in addition to the higher resolution setups 
used to achieve the blue and red parts of the stellar spectra. 
Careful attention was also paid to the cleaning of telluric ab¬ 
sorptions that were present in the redder part of the spectra. 


the original library since we removed those stars whose spec¬ 
tra were not properly representative of a given set of atmo¬ 
spheric parameters. These included spectroscopic binaries 
(identihed with SIMBAD), or stars with unexpectedly high 
signal of variability (a ccording to the CDS electronically- 
readable version of the iKholopov et al.lIl998l catalogue) for 
which there were alternative stars in the MILES database 
with similar parameters. Other reasons for removing stars 
were the presence of emission lines, very low S/N in the blue, 
problems in the continuum and the lack of atmospheric pa¬ 
rameters determinations. For this process, each stellar spec- 


Our models make use of a subsample of 925 stars from 
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Figure 2. Comparison of the PadovaOO and BaSTI isochrones for three representative metallicities (Z= 0.001, Z= 0.004 and solar), as 
quoted within the upper panels, and three different ages, 10, 1 and 0.1 Gyr from top to bottom panels. Note that the reference solar 
metallicity value is slightly different for these two models. There are differences in the RGB, HB and AGB phases as well as in the very 
low-mass dwarf stars. 


trum was compared to a synthetic spectrum of similar at¬ 
mospheric parameters computed with the interpolation al¬ 
gorithm described in Paper I after exclusion of the targeted 
star. Finally, 75 stellar spectra with some, milder, problems 
were retained in the final sample as they were found to be 
useful for improving the coverage of certain regions of the 
parameter space. For these stars we decreased their relative 
weight in comparison to the other stars. 

The parameters for MILES stars were obtained from 
an extensive compilation from the literature, homogenized 


by taking as a refere nce the stars in common with the 
high resolution stars o f ISoubiran. Katz fc Cavrell d 19981 ) (see 
ICenarro et al.l l2007ah . One of the most relevant features 
of this library is the significant improvement of the pa¬ 
rameter coverage in comparison to other widely employed 
stellar spectral libraries (see figure 2 in Paper I). MILES 
spans a fairly wide range in metallicity, including stars with 
[M/H]< —2 as well as super-solar, for the temperature range 
that is relevant for intermediate and old stellar populations, 
both for dwarf and giant stars. Note that, unlike in other 
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stellar libraries (see Paper I), MILES includes a number of 
metal-poor dwarf stars with Tefj> 6000 K, which are rel¬ 
evant for predicting stellar populations with ages smaller 
than 10 Gyr at low metallicities. This library also includes 
a number of very low mass dwarfs, which are relevant when 
varying the IMF slope. Although not optimised for very 
young stellar populations (see Paper I) the library has MS 
stars as hot as ~ 30000 K, mostly around solar metallicity. It 
also includes giants with temperatures above 6000 K, such as 
metal-poor Horizontal Branch (HB) stars, Post AGB stars 
and variable stars of the type <5 Scuti, RR Lyrae or Cepheids 
along the instability strip. 

An important improvement here is the consideration of 
the [Mg/Fe] estimates for most of th e MILES star s, whic h 
have been recently determined by iMilone et al.l ll201lll . 
These abundance ratios were obtained through a com¬ 
pilation of values from the literature using abundances 
from high-resolution spectroscopic studies and from a ro¬ 
bust spectroscopic analysis using the MILES spectra. These 
[Mg/Fe] values were carefully calibrated to a single uni¬ 
form scale with a typical uncertainty on the estimates of 
~ 0.1 dex. As expected, the MILES library foll ows the abun¬ 
danc e pattern of the Galaxy (see figure 10 in iMilone et al.l 
[ 2 oT 3 ), although there is some dispersion at any given metal¬ 
licity. There were 75 stars in th e selected MILES subsam¬ 
ple for which iMilone et al.l (l201lh could not determine their 
[Mg/Fe] values (the majority of them have effective temper¬ 
atures either above 8000 K or below 3500 K, and also 
stars belonging to stellar clusters). To retain these stars for 
our model computations we assigned the mean abundance 
ratio expected for th e [Fe/H] value of th e star according to 
the pattern found by IMilone et al.l || 2011 ^ for the field stars. 
Fo r the cluster s t ars w e adopted the mean value determined 
by IMilone et al.l ll 201 lli for the other stars of the same clus¬ 
ter. 

Fig. [3] shows the stellar parameter coverage of the sub¬ 
sample composed of 925 stars selected from the MILES 
database. Stars with [Mg/Fe] abundance ratio smaller 
(larger) than -1-0.2 are plotted in green (red). Note that 
we only plot parameters for stars with temperatures be¬ 
low 8000 K, as the abundance estimates for higher tem¬ 
perat ures is more uncertain (see for details IMilone et al.l 
[ 20 i 2 ) and the net effec t of the abundance on the spectra 
is small (see figure 9 of ISansom et al.ll2013li . As expected, 
the stars enhanced in [Mg/Fe] also have [Fe/H]< —0.5, 
whereas the higher metallicity regime is populated with 
scaled-solar stars. We also see that MILES stars show a 
good coverage of temperature within the metallicity regime 
—0.5 ^ [Fe/H] ^ -1-0.2. Despite a significant decrease in 
the number of stars for the more metal-poor and metal-rich 
regimes, there are sufficient stars to compute models with 
[Fe/H]~ —2.3 and [Fe/H]~ -1-0.4, although within a lim¬ 
ited age range. A quantitative analysis of the quality of the 
models based on the atmospheric parameters coverage of the 
library stars is shown in Section [4. II 


2.3.2 Theoretical stellar spectral library 

Theoretical spectra are employed to determine a differen¬ 
tial correction to be applied to the reference models com¬ 
puted on the basis of MILES (see Section |3}. In order 
to compute these “corrected” SSPs we adopt the theo- 



[Fe/H] [Fe/H] 


Figure 3. Stellar parameters coverage of the MILES library. Red 
triangles and green diamonds represent stars with [Mg/Fe] values 
larger and smaller than 0.2, respectively. For this plot we split 
the database in giants (left panel) and dwarfs (right panel) at 
log 3 = 3.0. 


retical stellar library of iGoelho et al.l (l2005l L with the ex¬ 
tension to very cool gi ants and the spect rophotometric 
calibration presented in iGoelho et al.l (l2007l L These stars 
were computed with Pfant spectral synthesis code, and 
are based on refined atomic and mol ecular line lists (e.g. 
ICavrel e t al.lll99 ll: iBarbuv et al.l[2003ll . Model atmospheres 
from Castelli &: KurucJ i 2003h are a dopted for ef fective tem¬ 
peratures Teff^ 3500 K, and from Pled ( |l99dl for Tefr< 
3500 K. We refer the reader to ICoelho et al.l 1 20051 . l2()07fl 
for more details. _ _ 

The IGoelho et al.l ll2005h theoretical library covers ef¬ 
fective temperatures, TeU, between 3500 and 7000 K, sur¬ 
face gravities, log 3 , between 4-0.0 and 4-5.0, iron metallic¬ 
ities, [Fe/H], between —2.5 and 4-0.5 and a-element over 
iron abundances, [a/Fe]= 4 -0.0 and [a/F e [= 4- 0.4. The ex¬ 
tension to cool giants from ICoelho et ^ (l2007l l covers Teu 
down to 2800 K, loggr between —0.5 and 1.5, [Fe/H]= —0.5, 
0.0 and 4-0.2, and [a/Fe]= 4-0.0 and [a/Fe]= 4-0.4. The 
scale d-solar spectra ([a/F e ]= 0.0 ) have the solar abundances 
from iGrevesse fc Sauval Ill998h . The a-enhanced spectra 
([a/Fe]= 4-0.4) adopt a uniform enhancement [X/Fe]= 4-0.4 
for the elements O, Ne, Mg, Si, S, Ca and Ti, and the other 
elements follow the solar abundances. 

The accuracy of t his library has been assessed in 
iMartins fc Goelhq (l2007l L who showed that this library had 
the highest success in reproducing spectral features of ob¬ 
served stars among the theoretical libraries analysed in that 
work. The spectra have higher spectral resolution and wave¬ 
length coverage than MILES, therefore they were convolved 
with a gaussian kernel and were resampled to match MILES 
spectral resolution, dispersion and wavelength coverage. 
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3 SSP SPECTRAL SYNTHESIS MODELLING 

To compute the spectra of single-burst stellar populations, 
i.e. SSPs, we follow the method described in detail in Pa¬ 
per I. This approach consists of integrating the spectra of the 
stars along the isochrone taking into account their number 
per mass bin according to the adopted IMF. In the present 
study we distinguish and provide three types of models for 
a given total metallicity, [M/H] 0. The first type of model 
SSPs are what we call base models, which combine scaled- 
solar isochrones of a given total metallicity with the MILES 
stars of similar [Fe/H] (rather than [M/H]), irrespective of 
their [Mg/Fe] abundance ratio. As MILES follows the abun¬ 
dance pattern of The Galaxy, these models are effectively 
scaled-solar at solar metallicity. However, for lower metallic- 
ities these models are not self-consistent in the sense that ob¬ 
served stars are enhanced in [Mg/Fe], whereas the isochrones 
are scaled-solar. Therefore the SSP SEDs of Paper I, which 
combine MILES with the PadovaOO isochrones, should be re¬ 
garded as base models. The second type of models presented 
here are termed scaled-solar as both the stellar spectra and 
the isochrones are scaled-solar throughout the whole metal¬ 
licity range. Finally, the third type of models computed here 
are the a-enhanced models, which combine a-enhanced stel¬ 
lar spectra with a-enhanced isochrones. For the scaled-solar 
and a-enhanced models, which are computed on the basis of 
the total metallicity, [M/H], we t ake into account bot h the 
recent ]Mg/Fe] determinations of iMilone et al.l ((2oT3) and 
the [Fe/H] estimates of the individual library stars. 


3.1 Base models 

The term ’’base model” refers to the commonly used ap¬ 
proach that is employed to compute SSP models, as we de¬ 
scribe in this Section. It should not be confused with the 
terms BaSTI- or PadovaOO-based models, which refer to the 
use of a particular set of isochrones used to generate the SSP 
models. For the base models we always employ scaled-solar 
isochrones. We compute here two sets of models that, to 
be precise, should be termed BaSTI-based base models and 
PadovaOO-based base models. For brevity henceforth we will 
refer to these models as BasTI base models and PadovaOO 
base models. 

The base SSP model spectra are computed according 
to the following expression 


]M/H] = ]Fe/H]. This is why we write in the left term of 
this equation [M/H], whereas in the right one we use the 
]Fe/H] abundance measurement. Therefore, the base mod¬ 
els in the low metallicity regime are not self-consistent since 
the MILES stars are enhanced in ]Mg/Fe], and therefore the 
true [M/H] value is in reality larger than the adopted [Fe/H] 
value (see Section[3]2|. A* (m, t) is the number of such stars 
in the mass bin (m, m-\-dm), which depends on the adopted 
IMF, <I>(m). The lowest and highest mass stars that are alive 
at this age are m\ and mt. Note that mt is determined by 
the isochrone. 

To normalise the stellar spect ra in the V band we con - 
volve with the filter response of iBuser fc Kuruc j lll978ll . 
These spectra are then scaled according to the predicted 
flux in this band, following the empirical photometric li¬ 
braries described in Section [221 In order to obtain the ab- 
solute flux in the E-band w e follow the method described in 
iFalcon-B arroso et al. ll2 011^. which is based on the calibra- 
tion of iFukugita. Shimasaku fc Ichikawal (Il995l i: 


_ 2^Q-o.4(v-i-Zpv-'4,i:,yr) 


( 2 ) 


where Zpv and VaLyr are the adopted zero-point and V 
magnitude for the reference Vega spe ctrum, respec tively. 
To compute the zero-point, we used the lHav^ ll 19851 ') spec¬ 
trum of Vega with a flux of 3.44 x 10“® erg cm“^ s“^ at 
5556 A. The V magnitude of Vega is s et to 0.03 mag , whic h 
is consistent with the value adopted bv lAlonso et al.l l|l995h . 
i.e. in agreement with our prescriptions to transform the the¬ 
oretical parameters of the isochrones to the observational 
plane. This scheme allowed us to achieve absolute magni¬ 
tudes derived directly from the synthesised SSP spectra, 
which were fully consistent with those computed from the 
photometric libraries for the same SSPs. These uncertainties 
are well below the observational uncertainties 0.01 mag) 
for all ages, metallicities, a nd IMFs. We refer the reader to 
iFalcon-Barroso et all (l201l[ ) for more details on the achieved 
accuracy. 


3.2 Truly scaled-solar and a-enhanced models 

To compute both the truly scaled-solar and the a-enhanced 
SSP spectra, with [a/Fe] = 0.0 and -1-0.4, respectively, we 
feed the code with stars according to their estimated [Fe/H] 
and [Mg/Fe] values. We therefore re-write Eq. [T] as follows 


SA(t,[M/H],-I>,7o.o) 


SAv.(m,t,]Fe/H]) x 


iV(m, t, [Fe/H])A#(m, t)dm, (1) 


where (”t-, t, ]Fe/H]) is the empirical stellar spectrum, 
irrespective of its [Mg/Fe] abundance, normalised in the 
V band, corresponding to a star of mass m and metal¬ 
licity ]Fe/H], which is alive at the age t of the stellar 
population. Its spectrum is also characterised by the Tefj 
and log g parameters as determined by the adopted scaled- 
solar isochrones, 7o.o- t, [Fe/H]) is the flux of the 

star in the V-band. For these base models we assume that 


® We use [M/H], rather than [Z/H] (as used in Paper I), because 
Z denotes a mass fraction, and therefore [M/H] is less confusing. 


pmt 

Sa(L [M/H], [a/Fe],-I>,7„) = / 5 a( m, t, ]Fe/H], [Mg/Fe]) 

J m\ 

x7^v(m,t, [Fe/H], [Mg/Fe]) 
xN,s,{m,t)dm, (3) 

In this case the isochrone, la, is characterised by the [a/Fe] 
abundance ratio, which in the present study can take the 
value of 0.0 and -1-0.4. This notation can be generalised to 
any abundance ratio, including any particular choice of el¬ 
ement ratios, by substituting a and ]Mg/Fe] by the desired 
elemental abundance(s) ratio(s) in this equation. Conversely, 
this notation is not appropriate for the base models, which 
do not take into account any consideration for the abun¬ 
dance ratios. 

Note that the only abundance ratio provided in 
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iMilone et alJ (l201lf) is that of [Mg/Fe], but not those cor¬ 
responding to the other a-elements. Hereafter we write 
[Mg/Fe] when referring to the abundance ratio of the stars, 
which is a measured quantity, whilst for the models we use 
[a/Fe]. This is because we use the [Mg/Fe] as a proxy for the 
[a/Fe] abundance ratio of the star (see below and see also 
Section 13.311 and because we employ theoretical isochrones 
with varying abundance for all the a-elements, i.e., [a/Fe]. 

We take into account the following relation between the 
iron abundance and the total metallicity 


[Fe/H] = [M/H] - A X [Mg/Fe] 


(4) 


The conversion between [Fe/H[ and [Mg/Fe] is dependent on 
the exact solar chemical mixture adopted. For the present 
study we adopt A = 0.75, in order to be consistent with the 
theoretical stell ar spectra of Section [T3^ which follow the 
solar mixture of iGrevesse &: Sauvall (Il998ll . 

For the high metallicity regime, i.e. with [M/H] 

0, where most of the MILES stars have [Mg/Fe] ~ 
0.0, we use Eq. [3] to combine MILES with scaled-solar 
isochrones to compute scaled-solar SSP spectra of total 
metallicity [M/H] = [Fe/H]. Therefore the resulting model, 
SA(t, [M/H], 0.0, <0, Jq.o), does not differ significantly from 
the corresponding base model, S\{t, [M/H], <I>), lo.o), around 
solar metallicity. Note that in this case only MILES library 
star spectra are employed to compute these models. 

To synthesise an a-enhanced SSP spectrum 
with [Mg/Fe] = -1-0.4 and total metallicity, [M/H], 

we first compute a reference SSP spectrum, 
[^/H]) O-Oi 4?,/ 0 . 4 ), which combines an a- 
enhanced isochrone, Jo. 4 , with MILES scaled-solar stars. 
We use the BaSTI database for this purpose. Note that 
this reference model is not self-consistent as the abundance 
ratios of the isochrones and the stellar spectra do differ. 
In addition, using the same a-enhanced isochrone we 
employ the theoretical stellar spectral library to com¬ 
pute two fully-theoretical SSP spectra: one scaled-solar 
and the other a-enhanced, (t, [M/H], 0.0,'F,/ 0 . 4 ) 

and (t, [M/H], 0.4, / 0 . 4 ), respectively. With these 

models we obtain a differential correction as follows 


7A(t,[M/H],0) 


gA..eo(^.[M/H],0.4,1>,7o.4) 

[M/H],0.0,-I>,/o.4) 


(5) 


This correction is then applied to the reference model, 
'S'xMri,Es(^> [WH],0.0,<E>,7 o.4), thereby giving a semi- 
empirical, MILES-based, self-consistenll^°l , a-enhanced SSP 
Spectrum 


It is worth noticing that, being similar, the theoretical stel¬ 
lar spectra adopt a uniform enhancement for all the a-elements, 
whereas the isochrones employ a different mixture for these ele¬ 
ments, which on average yields the same enhancement. The vari¬ 
ation of the abundance of individual a-ele ments has been shown 
to affect the temperature of the stars fsee lLee et al. I 120101 for an 
extended discussion). In addition, as described in Section 12.3.11 
we only take into consideration the [Fe/H] and [Mg/Fe] (which is 
used as a proxy for the overall a-enhancement) estimates that are 
available for the MILES stars. Therefore our models cannot be 
regarded as fully self-consistent as this is only possible for fully 
theoretical models that combine theoretical stellar spectra and 
isochrones computed with exactly the same element partition. 
With these important caveats in mind, our models are hereafter 
regarded as self-consistent (although not fully self-consistent) in 



Figure 4. The method followed for computing an a-enhanced 
SSP spectrum with [M/H] = -h0.06. Fully-theoretical a-enhanced 
([Mg/Fe] = -hO.4; shown in red) and scaled-solar (shown in green) 
SSP spectra, both computed with an a-enhanced isochrone, are 
employed to obtain a differential correction (in blue). This cor¬ 
rection is applied to a reference SSP spectrum (also shown in 
blue), computed with MILES scaled-solar stars and the same a- 
enhanced isochrone. Therefore we obtain a self-consistent SSP 
spectrum (shown in black) that is based on MILES, yet incor¬ 
porates the [a/Fe] abundance ratio effects through this theoret¬ 
ical differential correction. All the models have [M/H] = -1-0.06, 
10 Gyr and Kroupa Universal IMF. See the text for a detailed 
description. 


SA(i, [M/H], 0.4, $,70.4) = SA„,^^s(t,[M/H],0.0,$,7o.4) 

X7A(t,[M/H],$). (6) 

The key steps of the method are illustrated in Fig. [4] 

For low metallicities, i.e. with [M/H] < 0.4, where the 
MILES stars have [Mg/Fe] ~ -1-0.4, we use Eq. [3] to com¬ 
pute fully-consistent models employing the MILES stellar 
spectra and a-enhanced isochrones, with no aid from the¬ 
oretical stellar spectra. The MILES stars are selected to 
have [Mg/Fe] = +0.4 and [Fe/H] = [M/H] - 0.3, accord¬ 
ing to Eq. 131 To compute the scaled-solar models we follow 
a similar approach to that employed above for obtaining an 
a-enhanced SSP spectrum at high metallicity. However, in 
this case we apply a differential correction, 7 a( 7, [M/H], $), 
that allows us to transform the reference MILES-based 


the sense that the average o-enhancement value of the main in¬ 
gredients of these models is similar. 
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model with [a/Fe] = +0.4, (^> 0.4, lo.o) 

to a self-consistent SSP spectrum with [a/Fe] = 0.0, 
SA(t,[M/H],0.0,-F,7o.o). 

For intermediate metallicities, namely [M/H] = —0.35 
and [M/H] = —0.25, where most of the MILES stars have 
[Mg/Fe] ~ 0.2, we obtain for each [M/H] metallicity and 
each [a/Fe] abundance two SSP spectra using the two ap¬ 
proaches described above. On one hand, we compute self- 
consistent scaled-solar SSP models that are fully based on 
MILES making use of the stars that have [Mg/Fe] ~ 0.0 and 
[Fe/H] = —0.35 and [Fe/H] = —0.25, respectively. Similarly, 
for the corresponding a-enhanced models we select stars 
with [Fe/H] = —0.65 and [Fe/H]= —0.55, respectively. No¬ 
tice that at such low [Fe/H] abundances we take the advan¬ 
tage that most of the stars of MILES have [Mg/Fe] ~ +0.4. 
On the other hand, for each of these two [M/H] and two 
[a/Fe] values, we also compute a reference model based on 
MILES as well as two models based on the theoretical stellar 
spectral library. The latter are used to obtain the differen¬ 
tial correction that allows us to bring the reference model 
onto the desired [a/Fe] abundance ratio. The use of these 
two self-consistent SSP models for each [M/H], [a/Fe] com¬ 
bination, allows us to significantly increase the total num¬ 
ber of the MILES stars that are employed for the compu¬ 
tation of the final SSP spectrum, which, therefore, has a 
higher quality than either of the above two SSP spectra, 
separately. It is worth reiterating that we employ this ap¬ 
proach as our isochrones only have two possible [a/Fe] val¬ 
ues, i.e. 0.0 and +0.4, whereas the MILES stars mostly have 
[Mg/Fe] ~ +0.2. Finally, the last step consists of combin¬ 
ing these two SSP spectra using equal weights. This equal¬ 
weighting is motivated by the fact that the quality parame¬ 
ters of these two SSP spectra (described in Section [4.HI are 
very similar for the vast majority of the cases, with differ¬ 
ences being smaller than ~ 10 %. 

Our approach differs from other methods that apply 
a differential correction on the individual stellar spectra 
befor e synthesising the SSP s pectra (e.g., IPrueniel et al.l 
I 2 OO 7 I : I Conroy ^ van DokkurnI l2012al L In our approach 
7 a( t, [M/H],</) is obtained from the SSP spectra calculated 
specifically for this purpose, without modifying the stel¬ 
lar s£ectraHri_oimn:iod^, simi lar to the method pr oposed 
in ICervantes et al.l (120071 ) and IWalcher et al.l ll2009l). This 
differs from the differential models built in Walcher et al.l 
1 I 2 OO 9 I) in that all our intermediate SSP spectra, either fully- 
theoretical or reference MILES-based, are computed with 
the same population synth e sis cod e. We also differ from 
e.g. Conroy fc van DokkumI ll2012ah and ICervantes et al.l 
ll2007l) . in that we employ a-enhanced isochrones to obtain 
the a-enhanced SSP spectra. Our m odels also diff e r from 
the fully-theoretical SSP SEDs of, e.g. ICoelho et all (l2007h . 
in that we always use MILES as our input library either 
completely, or in the form of the reference models. To our 
knowledge, this is the first time that SSPs using the differ¬ 
ential approach and empirical stellar spectra have been built 
self-consistently. 

Fig. O shows the flux-ratios obtained by dividing fully- 
theoretical (i.e., no MILES stars have been used here) SSP 
spectra versus the corresponding models computed with 
MILES (i.e., with no aid of theoretical stellar spectra) for the 
same SSP parameters. These ratios correspond to two rep¬ 
resentative metallicity and abundance ratio regimes, where 



4000 5000 6000 7000 

Wavelength (A) 


Figure 5. Flux-ratio obtained by dividing an SSP SED computed 
with theoretical stellar spectra by the corresponding SSP model 
computed with MILES stars. In the top panel the models are 
scaled-solar and have [M/H] = +0.06, whereas in the lower panel 
the models are a-enhanced with [M/H] = —0.96. All the models 
have 10 Gyr and a Kroupa Universal IMF. 


it is possible to compute self-consistent SSP models based 
solely on the MILES stars: low-metallicity, a-enhanced and 
high-metallicity, scaled-solar. The two flux-ratios show a 
large-scale variation of the continuum in the blue spectral 
region. These residuals include relevant features, which are 
related to carbon and some also to nitrogen, such as CNO at 
~ 3860 A, Call H-K lines around ~ 3950 A, CN at 4150 A, 
Ca at ~ 4227 A, CH at ~ 4300 A and C 2 at ~ 4670 A (see 
iTripicco fc Belllll995h . Note that all these residuals are em¬ 
phasised in the flux-ratio obtained for the high metallicity 
regime. These features are attributed, at least in part, to 
possible differences in the CN abundance between the the¬ 
oretical and the MILES stars. Note that the CNO group is 
an important contributor to the total metal content and in¬ 
tegrated opacity in a stellar photosphere, whose patterns do 
not necessarily follow the a-elements, iron-peak elements or 
global me tallicity, but have thei r own significant contribu - 
tions ('e.g. lMcWilliam et al.ll2008h . [Martins k. Coelhol (l2007f) 
find for the C-dominated features that the difference be¬ 
tween the theoretical and empirical stars is emphasised for 
stars cooler than ~ 4500 K. These stars have a large impact 
on the old SSPs represented in Fig. [5] 

We also see in Fig. [5] significant residuals for all the 
hydrogen lines, e.g. H5 at ~ 4100 A, H 7 at ~ 4340 A, H/I at 
~ 4860 A, Ha at - 6563 A. This is a well known problem 
in the theoretical models; the modelling of these lines with 
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the assumption of LTE is able to reproduce the wings, but 
cannot match the core of the Balmer lines which form in the 
chromosphere and are not accounted for in the theoretical 
stellar spectra fsee lMartins &: CoelhollioOTl ). 

Residuals for many metallic lines are seen along the 
whole spectral range in the high metallicity regime. There 
are also large-scale residuals red-ward of 6000 A, which 
are associated with molecular bands (mainly TiO). The con¬ 
tribution from these bands is only relevant for cool stars 
4500 K), but is a contribution that increases very rapidly 
with decreasing temperature. The modelling of these bands 
is challenging in part due to an incomplete knowledge of 
the molecular (and line) lists. Difficulties also arise since 
there may be differences between the temperature scales 
of the theoretical and e mpi rical stars (s ee a discussion in 
iMartins &: Coelh3l2007l an d[c oe ih3l2014l) . 

We are aware that our method is not completely free 
from these residuals that are related to the use of the theo¬ 
retical stellar spectra. However, their effects on our models 
are minimised as we only employ a differential correction 
in passing from, e.g., a scaled-solar to an a-enhanced abun¬ 
dance ratio. This approach avoids most of the uncertain¬ 
ties coming from the prescriptions adopted in the theoreti¬ 
cal models such as the treatment of convectio n, turbulence, 
N-LTE effects, incomplete line lists, etc. See ICoelho et al.l 
(1200711 for a detailed discussion on these differences, includ¬ 
ing a comparison of the resulting line-strength indices. 

Finally, the very large residuals seen at wavelengths 
smaller than 3800 A are artificial in the sense that they 
are a result o f the s tellar spectra calculations performed by 
ICoelho et al.l (l2005l l who only included the first ten Hydro¬ 
gen Balmer transitions. This results in artificially high fluxes 
in the wavelength range AA 3650 - 3790 A in stars of spectral 
types F and G, and artificial peaks centred around 3700 A 
in the differential spectra computed from integrated light 
models. Obviously this problem is far more severe for the 
young and/or metal-poor stellar populations. For this rea¬ 
son we do not recommend the use of the spectral range blue- 
ward of ~ 3800 A for our model predictions that make use of 
the differential correction, i.e. the scaled-solar SSPs at low 
metallicities, and the a-enhanced SSPs at high metallicities. 


3.3 Stellar spectra interpolation scheme 

To compute each requested star spectrum, following Eq. [T] 
we use the local i nterpolation algorithm described in 
IVazdekis et alj (l2003tl (Appendix B) and Paper I. For a given 
set of stellar atmospheric parameters, 9o (i.e. 5040/Teff), 
log go, [Fe/H]o, we interpolate the spectra of adjacent stars. 
The code identifies the MILES stars whose parameters are 
enclosed within a given box around the point in parame¬ 
ter space. If no star is found this box is enlarged in each of 
these directions until suitable stars are found. This is done 
by dividing the original box into eight cubes, all with one 
corner at that point. This scheme is particularly useful for 
minimising biases in presence of gaps and asymmetries in 
the distribution of stars around the point. The larger the 
density of stars around the requested point, the smaller the 
box. The sizes of the smallest boxes are determined by the 
typical uncer tainties in the determin ation of the parameters 
according to ICenarro et all ll2001bl L In each of the boxes, 
the stars are combined taking into account their parameters 


and the S/N of their spectra. Stars with similar parame¬ 
ters and higher S/N are given more weight according to a 
gaussian function that is applied on each parameter axis. 

Finally, the combined spectra in the different boxes, 
with stellar parameters close to the requested ones, are used 
to obtain the final MILES-based synthetic spectrum. In or¬ 
der to obtain the required atmospheric parameters we em¬ 
ploy a system of linear equations. Note that in some cases, 
usually in scarcely populated regions of parameter space, 
when a mathematical but unphysical (i.e. with at least one 
star weighting negatively) solution is found we simply com¬ 
bine these stars according to the same weighting scheme 
applied previously for the individual boxes. This generally 
provides a star with parameters v ery close to tho s e init ially 
requested. We refer the reader to IVazdekis et al.l (I2nn3l '> for 
a full description of the algorithm. In the present work we 
have further minimised these unphysical cases by repeating 
the whole process after varying the initial set of parameters 
{6o and log go) in steps of IK and 0.1 dex for temperature 
and gravity, respectively, until a solution is found. We only 
allow a very modest deviation of these two parameters, i.e. 
10% of the initial box size, as estimated for that parametric 
point. We do not vary [Fe/H]o as in doing so we could obtain 
solutions that would be biased towards the solar metallicity 
direction, particularly for the lowest and highest metallicity 
bins covered by the models. 

We have extended this 3-dimensional interpolation al¬ 
gorithm to be able to synthesise the star spectra taking into 
account the [Mg/Fe] parameter, as a proxy for [a/Fe]. How¬ 
ever we are constrained within the limits imposed by the 
library coverage of this parameter. For this reason, follow - 
ing the abundance pattern of MILES l|Milone et ahllioill L 
at low metallicities our a-enhanced models are fully based 
on MILES, whereas for obtaining truly scaled-solar mod¬ 
els we require first MILES-based reference models (that 
are then modified by theoretical differential corrections). At 
high metallicities we make use of the reference models to 
obtain our a-enhanced predictions, whereas the scaled-solar 
models are fully based on MILES. Therefore our approach, 
which makes use of the [Fe/H] and [Mg/Fe] estimates avail¬ 
able for the stars of MILES, is fully optimised to exploit the 
potential of this library (and other empirical libraries). 

To synthesise a stellar spectrum with given [Mg/Fe] o 
and [Fe/H]o, we run our 3-dimensional interpolator twice: 
in one case, a;, we restrict the database to stars with 
[Mg/Fe] ^ [Mg/Fe] ^-I-A [Mg/Fe] and in the other, ah, we use 
the stars with [Mg/Fe] ^ [Mg/Fe]^ — A[Mg/Fe]. A[Mg/Fe] 
represents a fraction of the typical uncertain ty in the de¬ 
termin ation of this parameter according to iMilone et al.l 
ll201lll . By adopting larger A [Mg/Fe] values we increase 
the number of stars involved in the process, but at the ex¬ 
pense of increasing the probability of obtaining two stars 
with lower or higher [Mg/Fe] values than requested. After 
performing several tests, where we varied A[Mg/Fe] for a 
representative set of stellar parametric regions, we adopted 
A[Mg/Fe] = -1-0.05. Note that this corresponds to approx¬ 
imately half the typic al uncertainty in t he determination 
of this abundance fsee lMilone et aLl201ll L To estimate the 
[Mg/Fe] abundance of these stars, i.e. a; and a/,, we take into 
account the [Mg/Fe] estimates for all the stars that were in¬ 
volved in their calculation through the 3-dimensional inter¬ 
polation scheme. We then apply the same weights obtained 
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during this process, which depend on their proximity to 6q, 
log go, [Fe/H]o. Note that these three stellar parameters are 
already in place for ai and ah through this interpolation. 
Finally we combine them to obtain the desired stellar spec¬ 
trum with [Mg/Fe]o. In those cases where ai and ah stars 
have either lower or higher abundance ratio than [Mg/Fe]o 
we choose the star whose abundance value is the closest 
to the requested one. We note that such cases only occur 
in poorly populated parametric regions fe.g. at metallicity 
[M/H]= +0A or [M/H]= -2.3). 


4 SSP SPECTRA 

Table [T] summarises the spectral properties of the newly syn¬ 
thesized SSP SEDs. Th e nominal resolution of the mod¬ 
els is FWHM= 2.51 A llFalcon-Barroso et al.l l201lh . This 
resolution is constant along the whole spectral range AA 
3540.5 — 7409.6 A. The models have flux-calibrated re- 
sponse and are cleaned from telluric absorption residuals 
dSanchez-Blazauez et al.ll2006l l. For computing the SEDs we 
have adopted a total initial mass of 1 Mq . See also the 
caveats in Table [T] and Section 13.21 regarding the spectral 
range blue-ward of ~ 3800 A. 

4.1 Quality of the SSP spectra 

Table [2] lists the covered SSP parameters and also gives the 
ranges where these predictions cannot be considered safe. 

To estimate the reliability of the synthesised SSP spec¬ 
tra we follow the quantitative approach described in Pa¬ 
per I. For each SSP calculation we use the same interpo¬ 
lation scheme described in Section E31 which synthesises a 
representative stellar spectrum for a given set of atmospheric 
parameters. We compute a normalised quality parameter Qn 
following Equations 4-6 of Paper I. Essentially, the higher 
the density of stars around the requested parametric point 
of each star along the isochrone, the higher the quality of 
the SSP spectrum. Qn is determined by the MILES library, 
either when the model is computed on the basis of this li¬ 
brary alone or when we make use of the MILES reference 
model. Our quality parameter is normalised with respect to 
a minimum acceptable value, which comes from a poor, but 
still acceptable, parameter coverage. We refer the interested 
reader to Paper I for a full description of the method. The 
only difference in the present work is that in the process of 
synthesising each requested stellar spectrum we obtain two 
values for this parameter, which correspond to the two stars 
required by our extended interpolation algorithm to match 
the [Mg/Ee] abundance ratio. We simply combine these val¬ 
ues in the same way we did for these stars. 

The reliability of the scaled-solar and a-enhanced SSPs 
with bimodal and unimodal IMFs is assessed in Fig. [6] and 
Fig.[71 respectively. Each panel shows how Q„ varies with the 
SSP age and total metallicity. SSPs with values above 
1 can be considered of sufficient quality, and therefore they 
can be safely used. We use the same scale for all the panels 
of these two figures. The figures show that the values 
corresponding to the scaled-solar and a-enhanced SSPs are 
very similar. The subtle differences seen between the upper 
and lower panels are attributed to differences in the temper¬ 
atures of the selected stars along the isochrones, i.e. scaled- 
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Figure 6. Reliability of the SSP spectra synthesised with a bi¬ 
modal IMF with slope F;, = 1.3 (i.e. Kroupa-like; left panels) 
and = 3.3 (right panels), for scaled-solar (upper panels) and 
a-enhanced (lower panels) models. The SSP age increases from 
left to right, whereas the total metallicity, [M/H], is shown as dif¬ 
ferent line-types and thicknesses as given in the first panel. Only 
Qn values above 1 are acceptable, as indicated by the horizontal 
black dashed-line. 


solar vs. a-enhanced), and the [Fe/H] metallicity according 
to Eq. m 

Although not shown, the Qn values of the SSP mod¬ 
els computed with a Kroupa Universal and Kroupa revised 
IMFs are similar to those obtained for a bimodal IMF with 
slope 1.3. The quality of the scaled-solar models is lower 
than that obtained for the base models, i.e. above 5 (see 
Figure 6 of Paper I). This is because in the computation of 
truly scaled-solar models we select stars according to both 
[Fe/H] and [Mg/Fe], whereas for the base models we only 
take into account [Fe/H]. In summary, the base models, ei¬ 
ther based on PadovaOO or BaSTI have in general higher 
quality than the corresponding scaled-solar and a-enhanced 
models. 

As expected, the higher Qn values are obtained at solar 
metallicity. Qn reaches a value of ~ 4 for stellar populations 
in the range 1-10 Gyr for the SSP models with standard 
IMF slope (i.e. the left panels of Fig. |S] and Fig. [7]). 

We also see that the quality of the bimodal models is 
slightly higher than that computed with a single power-law 
IMF. This is due to the small number of very low-mass 
dwarfs with temperatures below ~ 3500 K. This also ex¬ 
plains the decay in quality for ages above ~ 10 Gyr, as 
well as the drop seen for the models with a bottom-heavy 
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Table 1. Spectral properties of the model SEDs. 


Spectral range 

AA 354 O. 5 I - 7409.6 A 

Spectral resolution 

FWHM= 2.51 A (cr = 64 kms-i at 5000 A) 

Linear dispersion 

0.9 A/pix 

Continuum shape 

Flux-scaled 

Telluric absorption 

Cleaned 

Units 

F;^/L0A“^Mq“^, Lq = 3.826 X 10®®erg.s“^ 


^ Unsafe blue-ward of 3800 A for the scaled-solar models with [M/H] < 0.06 and o-enhanced models with [M/H] > —0.4 (see Section l3.2ll . 


Table 2. SSP parameter coverage. 


IMF 


Mass range 

IMF shape 

Unimodal IMF slope 
Bimodal IMF slope 

0.1 - 100 M0 

Unimodal, Bimodal, Kroupa Universal^, Kroupa Revised^ 

0.3 ^ F ^ 3.3^ (Salpeter IMF ^ if F = 1.3) 

0.3 ^ Ft, ^ 3.3 

Age range (Gyr) 

BaSTI-based models 
PadovaOO-based models 

0.030 - 14.0 

0.063 - 17.8 

Total metallicity ([M/H]) 

BaSTI-based models 
PadovaOO-based models 

-2.27"*,-1.79®’®’’^,-1.49®’®’i°’“’i^,-1.26i®’i"‘'i®’i®’i’^,-0.96i®,-0.66,-0.35,-0.25,+0.06,+0.15,+0.26,+0.40"^ 
-2.32"*,-1.711^’20,-1.31,-0.71,-0.41,+0.0,+0.22 

Abundance ratio ([a/Fe]) 

Base models 
scaled-solar^^ 
a-enhanced^ ^ 

MILES pattern 

+0.0 

+0.4 


^ Same unsafe models as for Bimodal IMF with F;, = 1.3 
^ Unimodal IMF models unsafe for F > 2.3 
^ Same unsafe models as for Unimodal IMF with F = 1.3 
^ Unsafe models 

^ Unimodal BaSTI base models unsafe for t < 0.09 

® Bimodal BaSTI base models unsafe for {t < 0.1,Ff, < 2.3), {t < 0.09,2.3 ^ F^, ^ 2.5), {t < 0.08,2.8 ^ F^, ^ 3.0), {t < 0.06,Fb = 3.3) 
^ scaled-solar and a-enhanced models unsafe 
^ Unimodal BaSTI base models unsafe for t < 0.07 

® Bimodal BaSTI base models unsafe for {t < O.OTjF^, < 2.8), {t < 0.06,F^, = 2.8), {t < 0.05,Ff, = 3.3) 

Bimodal scaled-solar models unsafe for {t < 0.08,F^j ^ 2.5), (t < 0.07,2.8 ^ F^, ^ 3.0), {t < 0.06,F(, = 3.3) 

Bimodal a-enhanced models unsafe for (t < 0.08,F(, ^ 2.8), (t < 0.07,3.0 ^ F^, ^ 3.3) 

Unimodal scaled-solar and a-enhanced models unsafe for t < 0.08 
Unimodal BaSTI base models unsafe for t < 0.04 
Bimodal BaSTI base models unsafe for (£ < 0.04,Fb < 2.8) 

Bimodal scaled-solar models unsafe for {t < 0.07,F^, ^ 1-8), (t < 0.06,2.0 ^ F^, ^ 3.0), {t < 0.05,F^, = 3.3) 

Bimodal a-enhanced models unsafe for (t < 0.07,F& ^ 2.0), (t < 0.06,2.3 ^ F& ^ 3.0), (t < 0.05,Fb = 3.3) 

Unimodal scaled-solar and a-enhanced models unsafe for t < 0.07 
Unimodal scaled-solar and a-enhanced models unsafe for t < 0.05 

Bimodal PadovaOO base models unsafe for (t < 0.075,F& < 2.5), {t < 0.07,2.5 ^ F^ ^ 3.0) 

Unimodal PadovaOO base models unsafe for t < 0.075 

See the caveats in Tabled for the spectral range blue-ward of 3800 A. 


IMF. Note that, as expected, this drop is more pronounced 
for the unimodal SSP models. Note also that the models 
with bimodal IMF and slope 3.3 yield very similar values 
to those obtained with unimodal IMF and slope 2.3. This 
is due to the fact that the relative fraction of dwarfs with 


masses below 0.5 M(7^ is very similar for these two IMFs (see 
iLa Barbera et al]l2013h . 

Finally the models with a unimodal IMF with F > 2.3 
are unsafe due to an emphasised contribution of very low- 
mass stars with temperatures cooler than 3500 K, which are 
scarce in the MILES library. We refer the interested reader 
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Figure 7. Same as in Fig. [6] for models with unimodal IMF with 
two slopes (F = 1.3, i.e. Salpeter, and F = 2.3). 


to lVazdekis et ahl (|2012| ) for the details regarding this lim¬ 
itation. Despite the claim of a steepening of the IMF slope 
with galaxy velocity dispersion, such extremely bottom- 
heavy IMF models do not seem to be required for fitting 
even the most massive galaxies, always wit h F < 2.3 (e.g. 
iLa Barbera et ah I I 2 OI 3 I : ISpiniello et al.l[20l3 ). 

The variation of the Qn values with age at a given 
metallicity originates from the varying number density of 
stars with temperature, such as the Turn-Off stars, which 
become cooler with increasing age. This is illustrated by 
looking at the number density of stars around solar metal¬ 
licity (i.e. along a vertical line in the right panel of Fig. 0. 
For ages smaller than 1 Gyr we also see a decrease in quality 
due to a decrease in the number of stars with high temper¬ 
atures. 

These figures also show that the derived quality de¬ 
pends on metallicity. In general, the more we depart from 
solar the lower the Qn value. Unlike the base models (see 
Figure 6 of Paper I) the SSPs with supersolar metallicity 
show higher Qn values than the models with [M/H]= —0.35 
and —0.25 shown here. This is because at these intermedi¬ 
ate metallicities the required reference models, either with 
[a/Fe]= +0.0 [M/H],0.0, <&,/o. 4 )) or [Mg/Fe] = 

+0.4 [^/H],0.4, <0,/o.o)), composed of sig¬ 

nificantly fewer stars, since the majority have [Mg/Fe] ~ 
+0.2. For the base model we neglect this parameter and 
therefore all these stars are taken into account in the inter¬ 
polation scheme. The drop in quality is more pronounced for 
lower metallicities. Overall, the very youngest models are not 


Figure 8. Flux-ratios obtained by dividing an o-enhanced versus 
a scaled-solar SSP model with constant iron content ([Fe/H] = 
—0.25, in blue), and constant total metallicity ([M/H] = +0.06, 
in black). Both flux-ratios correspond to models with 10 Gyr and 
Kroupa Universal IMF. 


reliable for metallicities below [M/H]~ —1.0. The minimum 
acceptable metallicity for the a-enhanced and scaled-solar 
models is reached for models with [M/H] = —1.49. Note 
that this limit is lower for the base models, which are safe 
for [M/H] = —1.79 (BaSTI) with ages above ~ 0.1 Gyr and 
for [M/H] = —1.71 (PadovaOO) with ages above ~ 0.08 Gyr. 
Finally we also see that the a-enhanced and scaled-solar 
predictions are unsafe for ages younger than ~ 0.07 Gyr. In 
Tabled we provide more detailed information on which re¬ 
gions of parameter space (age, metallicity, IMF slope and 
IMF shape) cannot be considered fully safe. 

Models outside the safe ranges discussed here should not 
be considered suitable for detailed spectroscopic studies that 
employ line-strengths or a full spectrum-fitting approach. 
However these models might be still useful, with caveats, 
for some applications such as, e.g., broad-band colours or 
for more complex populations, which are dominated by SSPs 
within the safe ranges. We also computed scaled-solar and 
a-enhanced models with [M/H]= +0.4 and [M/H]= —2.3, 
which are unsafe, but that could help us to visualise trends 
with metallicity of e.g. line-strengths. 
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Figure 9. Flux-ratios obtained by dividing two SSP spectra com¬ 
puted with the same isochrones but in one case employing a- 
enhanced spectra and in the other scaled-solar stellar spectra, 
shown in black. The ratios resulting from dividing two SSP spec¬ 
tra computed with the same stellar spectra but varying the abun¬ 
dance of the isochrones are shown in blue. The results for two total 
metallicity values (+0.06 and —1.26) and two ages (2 and 12 Gyr) 
are shown as quoted within the panels. 


4.2 Effects of varying the abundance on the SSP 
spectra 

The effects of varying the a-elements abundance ratio at 
constant total metallicity, [M/H], and at constant [Fe/H] 
abundance can be assessed from the flux-ratios shown in 
Fig. [8l When the total metallicity is held constant we see 
a significant increase in flux at wavelengths shorter than 
about 4500 A. This behaviour is also seen when dividing 
stellar spectra with varying [Mg/Fe], both for empirical stars 
from the MILES library and for theoretical stellar spectra. 
We refer the reader to the extensive study performed in 
Saj^om et al.l (l2013l ) for more details (see also lCassisi et all 
2004h . This excess of blue flux is attributed to the fact 
that the a-enhanced element mixture has a decreased iron 
abundance, and therefore lower opacity, with respect to the 
scaled-solar model of similar total metallicity (irrespective 
of the exact value of the term multiplying the [Mg/Fe] abun¬ 
dance in Eq. 0. This is confirmed by the modest excess in 
blue flux seen in the bottom spectral ratio of Fig. [8] which 
corresponds to models with similar [Fe/H] abundance. 

In addition, we see similar spectral features blue-ward 
of ~ 4500 A in both the spectral ratios, although gener¬ 
ally these are less pronounced when [Fe/H] is held constant. 



4000 5000 6000 7000 

Wavelength (A) 


Figure 10. Flux-ratios obtained by dividing a-enhanced versus 
scaled-solar SSP spectra of 2 Gyr (green) and 12 Gyr (black). All 
the flux-ratios correspond to models with solar (total) metallicity 
and Kroupa Universai IMF. 


These residuals are mainly related to the Call H-K and Cal 
features around ~ 3950 A and ^ 4227 A, respectively, the 
CN molecular absorptions at ~ 4150 A and in the complex 
blue-ward of ~ 3900 A and the CH at ~ 4300 A and also 
within this blue complex. These residuals likely originate in 
intrinsic differences between the C, N and O and Ca abun¬ 
dances bet ween the MILES sta rs and the theoretical stellar 
spectra of ICoelho et al.l |200^. In our SSP model compu¬ 
tation we only have considered the [Mg/Fe] and [Fe/H] es¬ 
timates for the MILES stars, as the measurements for the 
other elements are incomplete for these stars. However there 
are indications of such varying abundance patterns among 
the MILES and the theoretical stars. For example, applying 
the re lation between [C/Fe] and [Fe/H] of iTakeda fc Hondal 
ll2005l) . i.e. [C/Fe]= -0.21[Fe/H]+0.014, derived for a sam¬ 
ple of 160 FGK dwarfs/subgiants, with metallicities around 
the solar value, we obtain [C/Fe]~ +0.07 for [M/H]= +0.06 
and [a/Fe]= +0.4, whilst the theoretical stellar spectra 
adopt a scaled-solar abundance for carbon. Moreover, un¬ 
like in the theoretical stellar spectra, for which [O/Fe] tracks 
[Mg/Fe], observationally the abundances of these two a- 
elements are slightly deco upled around solar metallicity (e.g. 
ISoubiran fc Girardli2005l ). In addition, the CNO group rep¬ 
resents an important fraction of the total metallicity, con¬ 
tributing significantly to the integrated opaci ty in a stellar 
photo sphere (see the extensive discussion in ISansom et al.l 
l2013h . 
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Figure 11. Flux-ratios obtained by dividing a-enhanced by 
scaled-solar SSP spectra with constant total metallicity. The flux- 
ratios are shown for four metallicity values: [M/H] = —1.79 
(black), —1.26 (blue),—0.35 (green) and -1-0.26 (red). All the flux- 
ratios correspond to models of 10 Gyr with a Kroupa Universal 
IMF. 

Fig. [8] also shows that the Fe lines at, e.g., ~ 5270 A 
and 5335 A, are deeper in the scaled-solar models than the 
a-enhanced models (positive residuals in the r atios plotted). 
This h a ppens also at f ixed [ Fe/H], as noted in iBarbuv et alJ 
(l2003tl ; Icoelho et alJ (l2005l ~) : Mg is an important electron 
donor and as such can act to lower the true continuum. 
Thi s effect is also seen for the Nal doublet a t ~ 5895 A (see, 
e.g.. lVazdekis et ^1 19971 : TCoelho et al.ll200'i^ l and, to a lesser 
degree, in the variations of the molecular features red-ward 
of ^ 6000 A. In general these features, although weaker, are 
also present in the flux-ratio obtained when [Fe/H] is kept 
constant. The exception here is the Mg feature at 5175 A, 
where the [Mg/H] abundance of the two SSP spectra differ 
by —0.3 dex. See Section [5] for more details. 

In Fig. [§]we separate the effects of the isochrones from 
those of the stellar spectra on our SEDs. We show the re¬ 
sults for two ages (2 and 12 Gyr) and two total metallicities 
with [M/H]= -1-0.06 and —1.26. For each age and metallicity 
we compute SSP spectra with two different isochrones (a- 
enhanced, and scaled-solar), but without changing the abun¬ 
dance ratios of the stellar spectra. In addition, we also show 
the flux-ratios obtained when using the same isochrones, 
but varying the abundance ratios of the stellar spectra. The 
flux-ratios obtained indicate that the effects of the stellar 
spectra are significantly more relevant than the isochrones, 
particularly for solar metallicity and old stellar populations, 


confirming the previous findings bv ICoelho et al.l (l2007l . Sec¬ 
tion 6.4 ) and I Lee et al.l (l2009l l. The impact of the stellar 
spectra is mostly seen in the excess of flux in the blue, the 
strengths of many features and, to a lesser degree, on some 
molecular bands. However the effects of the isochrones are 
mostly reflected in the continuum over the whole spectral 
range covered by the models. Although smaller than the ef¬ 
fects of the stellar spectra, the impact of the isochrones is 
significant for the SSPs of 2 Gyr. Note also the small differ¬ 
ence in flux between the two flux-ratios. 

Fig. HOI shows the flux-ratios obtained by dividing SSP 
spectra with the same total metallicity (solar), but with two 
different [a/Fe] abundances, for two ages (2 and 12 Gyr). We 
see that the impact on the spectra due to the [a/Fe] abun¬ 
dance is significantly greater at old ages. Such behaviour 
is expected since the SSPs with 2 Gyr include hotter stars, 
for which the effects of varying the [g/Fel abundance are 
considerably smaller (see Figure 9 of ISansom et al.ll2013l l. 

Finally, Fig. [11] illustrates the effects of varying the 
[a/Fe] abundance as a function of the total metallicity, 
[M/H], while keeping constant the age. We see both a larger 
excess in blue flux and a larger line feature residuals, for 
higher metallicities. As was the case when varying the SSP 
age, this behaviour is also expected when varying the metal¬ 
licity since the stellar temperatures decrease with increasing 
metallicity. 


5 LINE-STRENGTH INDICES 

In this section we study the behaviour of the line-strength 
indices predicted by our new models. We measure the in¬ 
dices on the synthesised SSP spectra at resolution 8.4 A 
(FWHM) to match the resolution of the LIS-8.4A system 
that was introduced in Pape r I. Unlike the Lick/IDS system 
dWorthev fc Ottaviani|[l997l l. all these indices are measured 
with the same resolution on a flux-scaled response. 

To illustrate the main behaviours and differences of 
our models we first compare a number of representative 
line-strengths from our base models computed with the 
PadovaOO and BaSTI isochrones. In the following subsec¬ 
tions we compare the predictions from our scaled-solar and 
a-enhanced models for a larger number of commonly used 
line indices, which include those of Lick. We have organ¬ 
ised these indices according to their sensitivities to relevant 
parameters of the stellar populations. 

5.1 Line-strengths of the base models 

Fig. [12] shows the predictions of some line-strength indices 
as a function of age for the base models, computed with 
BaSTI and PadovaOO isochrones. The predictions are shown 
for three different metallicities. As can be seen, the differ¬ 
ences found using these two sets of isochrones are very small. 
For old stellar populations the BaSTI models give index val¬ 
ues that are slightly higher for the Balmer lines (around 
0.1 —0.2 A for old stellar populations), and slightly lower for 
the three metallicity indices shown in the figure. Note that 
the given [M/H] values are slightly higher for BaSTI than 
for PadovaOO. Therefore when comparing the BaSTI model 
predictions to observational data we expect the derived ages 
and metallicities to be slightly higher than when using the 
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Figure 12. Base models with two different isochrones; black 
are BaSTI and magenta lines are PadovaOO. Three differ¬ 
ent metallicities are plotted for BaSTI (PadovaOO), [M/H] = 
—0.35(—0.40),-1-0.06(0.00), and -h0.26(-|-0.22), with dashed, solid 
and dotted lines, respectively. All the models have a Kroupa Uni¬ 
versal IMF. Note that whereas the models based on the BaSTI 
isochrones reach ages as young as 30 Myr the PadovaOO-based 
models stop at 63 Myr (see Table n. All the index measure¬ 
ments were obtained in the LIS-8.4A system (Paper I), i.e. at 
FWHM= 8.4 A and flux-calibrated response. 



•Age=14 Gyr 


■Age=7.5' Gyr' 


■Age= 2 Gyr tM/H]=-o.66' 

[M/H]=-0.25 

[M/H]=0.06 


0.2 0.4 


[a/Fe] 

Figure 14. [MgFe] (left panels) and [MgFe]' (right panels) indices 
as a function of [o/Fe], for SSP models with three different ages 
(14, 7.5 and 2 Gyr) from top to bottom panels, and four differ¬ 
ent metallicities ([M/H]= —0.66,—0.25,-^0.06 and -hO.26, in red, 
green, blue and cyan, respectively. All the models have a Kroupa 
Universal IMF. Note that the two indices are almost insensitive 
to this abundance ratio. 


PadovaOO model predictions. There are also differences for 
intermediate-aged stellar populations where the AGB con¬ 
tribution is relevant, due to the differences in the treatment 
of this late evolutionary stage (see Section E3I- 

5.2 Line-strengths with varying abundance ratios 

Fig. [T3 shows the predicted time evolution of our BaSTI- 
based SSP models for a selection of line-strength indices 
measured from the scaled-solar, a-enhanced and base mod¬ 
els. All these models are computed with a Kroupa Universal 
IMF. We have plotted the results for four different metallic¬ 
ities. Perhaps most readily seen in this set of plots is that, 
for the majority of indices, the effect of the different [a/Fe] 
ratios is largest at old ages. For some indices (e.g., Mgi, 
Fe5015, Fe5335, Fe5406, Fe5709, Fe5782), differences also 
appear at very young ages (< 0.1 Gyr), with these differences 
increasing at high metallicities. These differences come from 
the influence of supergiant stars at these young ages. 

The base models have the abundance pattern of the 
Milky Way (MILES) stars, although the isochrones are 
scaled-solar. Therefore, at low metallicities ([M/H] < —0.4), 
the base models are more similar to the a-enhanced models, 
while at high metallicities they are almost identical to the 
scaled-solar models. 

As expected those indices sensitive to Fe (Fe4383, 
Fe4531, Fe5015, Fe5270, Fe5335, Fe5406, Fe5709 and 
Fe 5782) show lower values in the a-enhanced models than 
in the scaled-solar models. This occurs because the mod¬ 
els are built at constant metallicity and, therefore, the a- 
enhanced models have a lower Fe content (i.e, [M/H]—0.3, 
according to Eq. 0. Note however that the effect on the 
Fe 5015 index is relatively smaller in comparison to the 


other Fe-sensitive indices, as a result of its sensitivity to 
Ti abundance , which is enhanced in the a -enhanced mod¬ 
els dLee et al.l l2009l : I Johansson et al.l I 2 OI 2 I I . The same ap¬ 
plies to Fe4531 for which we also see that the a-enhanced 
model lines are closer to the correspond ing base model 
lines (see also figure 9 in iLee et ahl l2009l i This beahviour 
is also seen for hxed [Fe/H], as Mg, i.e an i mportant elec¬ 
tron donor, tends to lowe r the true continuum (iBarbuv et al.l 
l2003l: ICoelho et al.ll2005h . 

It is worth noting that, while Mgb is higher in the a- 
enhanced models than the scaled-solar models, as expected 
due to its sensitivity to the Mg abundance, the same is not 
true for Mgi and Mg 2 . In fact, the Mgi indices are stronger 
in the scaled-solar models at all metallicities. It is known 
that Mgi has a greater sensitivity to C than Mg. In the a- 
enhanced models, all the non-a elements are depleted rela¬ 
tive to the scaled-solar models, including C. This behaviour 
of Mgi is similar to that of the CN 2 and C 2 4668 indices, 
which have higher sensitivities to variations of N and C 
abundances, respectively. 

Similar to the case for Mg6, Ca4227 is also much 
stronger in the a-enhanced models than in the scaled-solar 
models, as expected due to the sensitivity of this index to 
variations of Ca. Notice that Ca4227 is also anti-correlated 
with the CN feature that affects its the blue pseudocontin¬ 
uum, hence the Ca4227 index decreases in strength in the 
scaled-solar models. As with the Mg6 index, the predictions 
of the a-enhanced models are similar to those of the base 
models at low metallicities. This is as expected due to the 
high [a/Fe] of the base models at low metallicities (reflect¬ 
ing the intrinsic MILES abundance pattern). The Ca4455 
index, however, does not depend on the Ca abundance, and 
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Figure 13. Variations of a representative set of Lick line-strength indices with age for four different metallicities, [M/H]=—1.26 (dot- 
dashed line); [M/H]= —0.35 (dashed line) [M/H]= 0.0 (solid line) and [M/H]= +0.26 (dotted line). Colours represent three different 
models built with the BaSTI isochrones; Red: a-enhanced models, green: scaled-solar models and black: base models. As in Fig. 1121 all 
the index measurements are performed in the LIS-8.4A system, i.e. at 8.4 A (FWHM) and flux-calibrated response. 


behaves more like an Fe-index (e.g., iTripicco fc Belli 1 19951 : 
iLee et al]l2009h . 

Finally, in Fig.[T4]we make use of these new model pre¬ 
dictions to probe the sensitivity of the combined [MgFe] and 
[MgFel' indices to the g-en hancement. We confirm the re¬ 
sults of lXhomas et al.l l|200,3h . These authors employed mod¬ 
els that predict the strengths of a number of line indices with 
varying abundances by implementing response functions de¬ 
rived from theoretical spectra. We conclude that these two 
indices can be safely used to estimate the total metallicity, 
[M/H], of stellar populations with non-solar [Mg/Fe] abun¬ 
dance ratios. 


5.3 Age-sensitive line-strength indices 

Fig. IlSI shows the variation of different age-sensitive indices 
with age predicted by our models. The sensitivity of the 
H/3 index to a-enhancement is much less than that of the 
higher order Balmer lines (H5 and Hy). This is particularly 
true at high metallicities, where the base, scaled-solar and 
a-enhanced models have similar values. The dependence of 
H/3o on [a/Fe] 0.2 A at solar metallicity and old ages) is 
higher t han that of H/3 (almost negl igible), confirming the 
result of ICervantes fc VazdekisI ll2009l ). Also the dependence 


of CallH-K on [a/Fe] is very strong, especially for old ages. 
This dependence also increases with metallicity. Note that 
Ca is another Q-element and, therefore, in the adopted the¬ 
oretical stellar spectra is also enhanced in the a-enhanced 
models. Similar behaviour is seen for the D 4000 break. 
These variations are expected given the residuals obtained 
in Fig. [8] (see also Section [4.21) . Note also the increase in the 
strength of all the Balmer indices for the models with the 
oldest ages and the lowest metallicity, i.e. [M/H] = —1.26. 
This effect is due to the increased contribution of the Hor- 
izont al Branch stars at such old ages and low metallicities 
(e.g., iLee. Yoon &: 1^120001) 


Fig. Uniillustrates the effects of a-enhancement on vari¬ 
ous Balmer line index definitions. The figure shows the pseu- 
doconti nua and feature bandpass def initions for the four in¬ 
dices of lWorthev fc Ottavianil (Il997l). (i.e. H^a, H(5f , H+a , 
H+f), one of the Hy indices of Vazdekis fc Arimotd (Il999l) 
(Hy(T 2 oo, i.e. the one that is suitable for the resolution of the 
LIS-8.4 A system) and the sta ndard and opt i mised index def- 
initions for the H;5 fea ture of iTrager et al.l lll998l) (H/3) and 
ICervantes fc VazdekisI (l2009l) (H/3o), respectively. 

Fig. [reconfirms the claim of iThomas et al.l (l2003l) that 
the higher order Balmer lines are more sensitive to the a- 
enhancement than H/3o. However we see that such an effect 
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Figure 15. Predicted variation of various age-sensitive line-strength indices in the LIS-8.4 A system as a function of age for four different 
metallicities, [M/H]= -1.26,-0.35,-1-0.06, and +0.26 (dot-dashed, dashed, solid and dotted line, respectively). Red line: a-enhanced 
models, green: scaled-solar models and black: base models. All the models have a Kroupa Universal IMF. 


is mostly caused by the sensitivity of the surrounding pseu- 
docontinua of the H(5 and Hy lines to the abundance ratio, 
which is significa ntly lower for t he YLP index (see also fig¬ 
ures 17 and 18 in lLee et al.ll200^ . Essentially, the narrower 
the index definition the smaller the sensitivity of the index 
to the a-enhancement. For example, HyA has a larger depen¬ 
dence on the a-enhancement than Hy_F, while Hy(j 2 oo has 
the lowest sensitivity to this parameter as shown in Fig. 1151 

Fig. 1161 also shows two small, but significant, residuals 
(between ~ 4850 A and ~ 4900 A) seen in emission in this 
flux-ratio, affecting the H/3 and H/3o inde x definitions (se e 
also the spectra of figures 12, 13 and 14 in lLee et al. lEoS). 
One of these residuals (i.e., the bluest) affects the feature 
bandpass of the two indices, whereas the second (the red¬ 
dest) is placed within the red pseudo-continuum correspond¬ 
ing to the H/3 index. We see that the very small sensitivity 
of the H/J index to the [a/Fe] abundance ratio is due to 


the fact that the reddest residual compensates in part the 
contribution to the feature bandpass from the bluest resid¬ 
ual. However, the red pseudo-continuum of the H/?o index is 
located in a valley within this flux-ratio, which has a sim¬ 
ilar flux level to the blue pseudo-continuum. Therefore the 
H/9o index includes the full effect of the residual affecting 
the band-pass of this index. This residual corresponds to 
Fe-peak elements, including Chromium, that affect the red 
wing of the H/3 absorption. 


The higher sensitivity of Hdo to [a/Fe] was first 
shown in Cervan tes fc VazdekisI (l2009l ). also employing the 
ICoelho et al.l ll2005l) model stellar spectra. However these 
authors also found the opposite behaviour, i.e. a larger sen- 
siti vity of the Hd inde x to [a/Fe] than H/9o, when using 
the lMunari et all (l2005l ) stellar spectra to obtain the differ¬ 
ential correction. This disagreement possibly originates in 
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Figure 16. Flux-ratio obtained by dividing an a-enhanced SSP 
model by a scaled-solar SSP model, with age lOGyr, [M/H]= 
-1-0.06 and Kroupa Universal IMF, in the spectral region covered 
by H(5, Hy and H/3. The flux-ratio is shown at two resolutions 
FWHM= 2.51 A (i.e. the nominal resolution of the models) and 
8 . 4 A (i.e. the LIS-8.4A system), in thin green and thick black 
lines, respectively. The pseudocontinua and feature bandpasses 
of various index definitions for these features are shown in blue 
solid and dotted lines, respectively. 



0.1 1 10 0.1 1 10 


Age(Gyr) Age(Gyr) 

Figure 17. Variation of some IMF-sensitive indices with age of 
SSPs with a bimodal IMF with slope F^ = 1.3 (left panels) and 
F^ = 3.3 (right panels) for the three different models presented 
here; a-enhanced (red), scaled-solar (green) and base (black). The 
index measurements are shown for three different metallicities 
with different line-styles: [M/H]= —0.35 (dashed), [M/H]= 4-0.06 
(solid) and [M/H]= 4-0.26 (dotted). All the index measurements 
are in the LIS-8.4A system 


differences in the line lists adopted for computing these two 
sets of theoretical stellar spectra. 

5.4 IMF-sensitive line-strength indices 

Figjm shows the variat ion of four gravity-sensitive indices 
dha Barbera et al.ll2013l ) with age for the three models pre¬ 
sented here (i.e. base, scaled-solar and a-enhanced models). 
The figure shows the index measurements for SSPs with 
a bimodal IMF with slope Fj = 1.3, similar to Kroupa 
Universal, and F;, = 3.3, i.e. a bottom-heavy IMF. These 
values cover the range of bimodal IMF slopes required to 
match SDSS ET Gs with velocity dispers ions in the range 
100 - 300 kms“^ dLa Barbera et cil]|2013l '). The index mea¬ 
surements are shown for three different metallicities. The 
two TiO indices show the highest sensitivities to the IMF, 
whereas the NaD index shows the lowest sensitivity to this 
parameter. However, the NaD index does show the largest 
sensitivity to [a/Fe], whereas the Ti02 index is almost insen¬ 
sitive to this abundance ratio. These results are in agreement 
with recent works that employ either a differ ent set of model 
stellar spectra than the ones used here, e.g.. ISpiniello et al.l 
(|2014| . or fully empiric al evidence based on sta cked galaxy 
spectra from the SDSS llLa Barbera et al.| [20l^). 

The highest sensitivity to the IMF is seen for the mod¬ 
els with the oldest ages. For models with Fb = 1.3, we 
find a plateau, with only a modest increase with age, for 
ages between 2 and 10 Gyr for the TiOi and Ti02 in¬ 
dices and, to a lesser degree, for Mg 4780. For this age 


regime these indices mostly depend on the IMF, whereas 
for older populations th ere is also a dependence on age (see 
iLa Barbera et aLll2013li . Interestingly, the a-enhanced mod¬ 
els show a slightly larger dependence on metallicity in com¬ 
parison to the scaled-solar and base models. We also find 
that the Ti02 index, and the other IMF-sensitive indices, 
which show little sensitivity to the age for the standard IMF 
slope (i.e 1.3), become more sensitive to this parameter for 
larger IMF slopes, such as those shown in the right panels 
of Fig. [iTl 

Another interesting aspect is the high values of the TiOi 
and Ti02 indices for ages below 0.1 Gyr for the models with 
Fi = 1.3, index values which are similar to those obtained 
for ages around ~8 Gyr. This behaviour is attributed to the 
influence of supergiant stars at these young ages. Such an 
increase is likely to continue to even younger ages where 
the contribution of this evolutionary phase peaks at around 
lOMyr. Similar behaviour is seen for other IMF-sensitive 
features such as the Call tr iplet around ~ 8500 A and the 
CO band-head at 2.3/r fe.g. lMavvalll997h . It is not surpris¬ 
ing to see that there is no such increase at these young ages 
for the bottom-heavy IMF models. The steeper the IMF, the 
larger the fraction of mass locked in low-mass stars. In prin¬ 
ciple, it could be appealing to consider that the high TiO 
values observed in massive ETGs, whic h suggest a bottom- 
heavy IMF (e.g. iLa Barbera et ^ IEqH), may result in part 
by a contribution from such young stellar populations. How¬ 
ever such a contribution would have a significant impact on 
the measurements of other line indices in the optical range 
which is not seen. 
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Figure 18. Dependence of broad-band colours on the isochrones adopted. Colours are displayed as a function of age for four different 
metallicities: [M/H]= —1.26 (dot-dashed line), —0.35 (dashed line), +0.06 (solid line) and +0.26 (dotted line). The magenta lines indicate 
the models computed with the PadovaOO isochrones, whereas the black lines show the models based on BaSTI. The B—V colour (upper-left 
panel) is computed in the Vega magnitude system, whereas the SDSS colours u—g (upper-right), g—r (lower-left) and r—i (lower-right) 
are computed in the AB system. The residual colours are also shown in each panel. All the models have a Kroupa Universal IMF 


6 COLOURS 


In this Section we focus on the photometric predictions from 
our new set of models. To derive the synthetic colours, we 
adopt Vega magnitude s for the Johnson bands, B and V 
(|Buser fc Kuruca Il978l : see also iFalcon-Barroso et al.l 1201 ll 
for further details on our definition of Vega magnitudes) 
and the AB system for the SDSS filters: u, g, r and i. 


6.1 Synthetic u and i magnitudes 

The MILES spectral range (3540.5 — 7409.6 A) does not 
fully cover the wavelength range spanned by the u and the i 
SDSS filters, as the u pass-band starts at ~ 3005 A and the i 
band extends up to 8605 A. Therefore, in order to have reli¬ 
able synthetic magnitudes in these bands, we need to correct 
for the missing flux blue-ward and red-ward of the MILES 
spectral range. We do this by using a simi lar methodology 
to that adopted in iRicciardelli et al.l (l2012l ). which consists 
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Figure 19. Left-hand panels: comparison between colours in the SDSS bands of scaled-solar models (green lines) and the base models 
(black lines). Right-hand panels: comparison between Q-enhanced models (red lines) and base models (black lines). The colour variation 
is shown as a function of age for four different metallicities: [M/H]=—1.26 (dot-dashed line), [M/H]=—0.35 (dashed line), [M/H]= -1-0.06 
(solid line) and [M/H]= -hO.26 (dotted line). The residual colors are also shown in each panel. 
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of computing the magnitude in the part of the filter covered 
by the model and then correcting it for the missing flux. 

To obtain the missing flux blue-ward of the starting 
point of MILES, i.e. between 3005 A and 3540.5 A, we use 
new model SEDs that employ the NG SL st ellar library, as 
prepared in iKoleva fc Vazdel^ (120121 ') and IVazdekis et al.l 
ll2014l l. These models, which are computed with the same 
code used here, extend the spectral range down to 1680 A. 
Thus, for a given MILES SSP, the u magnitude is obtained 
from the flux through the u filter red-ward 3540.5 A and 
then correcting this value by the missing flux obtained from 
the corresponding NGSL-based SSP. Note that the missing 
flux correction in the blue is < 25 % for most of the SSPs 
covered here. Strictly speaking, this flux correction is only 
valid for the base models, which are of the same type as these 
newly synthesised NGSL-based SSP spectra. However the 
difference in flux due to abundance variations is far smaller 
than our correction for the missing flux. 

The same approach is used for calculating the miss¬ 
ing flux in the i filter, but this time using our pub - 
lished models extension, MIUSGAT (IVazdekis et al.|[2(il2ll . 
For the MIUSGAT models we employ the Indo-US stel¬ 
lar library IValdes et al.l (120041 ') to fill-in the gap left in 
betwee n the MILES range (Paper I) and Call triplet 
range d'Vazdekis et al.l l2003l L Since the effects of the a- 
enhancement on the continuum red-ward the 'V band is small 
(as shown by the flux-ratios of Fig. HOI and Fig. 1111) . the ex¬ 
tension based on MIUSGAT is appropriate for this specific 
purpose. 



Figure 20. Ratio of a-enhanced to scaled-solar models, 
at fixed [Fe/H] = 0 . 0, for our models (green) and the 

IConrov fc van DokkumI ll2012al) (red) SSP models with an age 
of 13.5 Gyr and a Salpeter IMF. The green curve has been arbi¬ 
trarily shifted downwards by 0.1 in order to allow a better visual 
comparison with the red curve. 


6.2 Base models with varying isochrones 

Fig. [TS] shows the comparison between the broad-band 
colours based on PadovaOO (magenta lines) and BaSTI 
(black lines) isochrones for the base models. Metal-poor 
models based on BaSTI tend to have bluer colours in all 
bands with respect to the PadovaOO models, by up to 0.2 mag 
for the lowest metallicity considered here ([M/H]= —1.26). 
For solar and supersolar metallicities, a significant offset be¬ 
tween the two sets of models is seen only in the r—i colour, 
where the BaSTI models are bluer by 0.02 mag for old 
populations and by ~ 0.04 mag for young populations. Over¬ 
all the agreement is good for the old stellar populations, 
whereas for younger ages, these colour differences reflect the 
differing prescriptions for AGB stars. 


6.3 scaled-solar versus a-enhanced SSP models 

The comparison between base, scaled-solar and a-enhanced 
models is shown in Fig. [T§] for the SDSS colours. No signif¬ 
icant differences are observed between the base models and 
the scaled-solar ones, with the exception of the u—g colour 
for the most metal-poor SSPs. Such difference must be at¬ 
tributed to the fact that at such low metallicity the MILES 
stars are a-enhanced. Gonversely, a-enhancement has a re¬ 
markable impact on the resulting broad-band colours, hav¬ 
ing a differential effect on the different spectral ranges (see 
also iGoelho et al.l l2007l) . The a-enhanced models are bluer 
in the blue (by ~ 0.2 mag in u—g and ~ 0.05 mag in g—r for 
old populations) and redder in the red (by ~ 0.02 mag in 
r—i) with respect to the base models, as one can also infer 


from the flux ratio of the two SSP spectra (Fig. 1111) . Note 
that the best agreement between the base and a-enhanced 
models is reached at low metallicities. 


7 COMPARISON WITH OTHER MODELS 


We compare the effect of varying [a/Fe] in our mod¬ 
els with predictions from other synthetic stellar pop¬ 
ulation models in the literature. We focus on some 
widely employed models t hat are based on empir ic al stel - 
lar spect ral libraries, i.e. iGonrov fc van DokkurnI (l2012ah 
[Gv pi2], [T ho mas. Maraston fc Johansson ( 201ll ) [TMJll] 
and ISch^onfl'^OOit) fS07l. Note that while CvD12 provide 
full SEDs at varying abundance ratios, TMJll and S07 pro¬ 
vide line-strengths of Lick spectral features. We note how¬ 
ever that whereas the predictions of TMJll are given at the 
nominal resolution of MILES, the models of S07 are pro - 
vided in the Lick/IDS system dWorthev fc Ottaviani|[r997l ). 

Fig. [20] compares the ratio of a-enhanced ([a/Fe] = 
+0.4) to scaled-solar SSPs for our models and those of 
GvD12, at a fixed age of 13.5 Gyr, and for the same 
(Salpeter) IMF. The GvD12 models are computed at fixed 
[Fe/Hj. Therefore, in order to perform a meaningful com¬ 
parison, we have computed our a-enhanced SSP for a total 
metallicity of [M/H] = +0.3 (i.e. [Fe/H] ~ 0.0, see Eq. |4|), 
and taken its ratio with the scaled-solar model. Note also 
that both our models and the GvD12 models have the same 
spectral resolution (i.e. that of MILES) in the spectral range 
considered here (from ~ 3500 to 7400 A). The overall effect 
of increasing a-enhancement is similar in the GvD12 models 
and our models, reflected by the fact that the flux ratios in 
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Figure 21. Line-strengths for different Fe lines (Fe 4383, Fe 5270, 
Fe5335, and Fe5015; from top to bottom, and left to right) 
are plotted as a function of Mg6 (A 5177 A). Black arrows, 
plotted with different line types, show the effect of changing 
[a/Fe], at fixed [Fe/H] = 0.0, for SSP models with the same 
age of 13.5 Gyr and the same (Salpeter) IMF , from different 
authors, i.e. this work (solid ), iThomas et al.l ll201ll. das hed), 
IConrov & van Dokkuiiil ll2012al . dotted*). and lSchiavorJ (l2007l . dot- 
dashed), as labelled in the top—left panel. The blue and magenta 
arrows show, as reference, the effect of decreasing age by 4 Gyr 
and increasing total metallicity by -1-0.3 dex (for [a/Fe] = 0.0), 
respectively, in our SSP models. 


the Figure have a similar large-scale behavior, and similar 
“absorption” and “emission” residual features (e.g., for Mg 
and Fe lines, respectively). While the average flux ratio is 
about one for CvD12 models, it is less than one for our mod¬ 
els. This occurs because we adopt a-enhanced isochrones for 
our a-enhanced model predictions, whereas CvD12 always 
use scaled-solar isochrones. Note however that the impact 
of the isochrones on the line indices is significantly less than 
the impact of the stellar spectra on the indices; see Fig. [9] 
It is evident that the ratios of different features in the 
two curves in Fig. [20] are not precisely the same, indicating 
that the two sets of models differ to some extent in the rel¬ 
ative effect of [a/Fe] on different spectral features. This is 
not surprising since the theoretical stellar spectra employed 
by CvD12 differ from those used in our models (see Sec- 
tion l2.3.^ . In particular, CvD12 obtained their differential 
corrections using theoretical stars computed with the AT- 
LAS12 code for model a tmospheres and the SYNTHE spec- 
tr um synthesis pa ckage (lKurucdll97fll . 1 199,31 : ported to Linux 
bv lSbordonel2003 ) and line list provided by Kurucz. Another 
relevant difference is that our differential corrections are ob¬ 
tained from SSP spectra computed with theoretical stellar 
spectra (see Section[3]2|, whereas the differential corrections 
of CvD12 were obtained from individual stars. Finally, as 
mentioned above, CvD12 employ scaled-solar isochrones for 
both their scaled-solar and a-enhanced predictions, whereas 
we employ a-enhanced isochrones for our a-enhanced model 


predictions. However, we note that the choice of isochrones 
has a small impact on the index predictions. 

Some of the model differences seen in Fig. are quan¬ 
tified in Fig. 1211 where we compare the effect of [a/Fe] on 
the Mg and Fe lines for our models and also for those of S07, 
TMJll and CvD12. Similarly, Fig. l22l shows the behaviour of 
the line-strengths of some representative Balmer lines (upper 
panels) and IMF-sensitive (lower panels) indices as a func¬ 
tion of the total metallicity indicator [MgFe]’, for the differ¬ 
ent models. Each arrow in a given panel represents an indi¬ 
vidual model, with the black arrows connecting scaled-solar 
to [a/Fe] = 4-0.4 SSP model line-strengths, at fixed ]Fe/H]. 
Since our models and those of TMJll are computed at fixed 
total metallicity, the comparison at fixed [Fe/H] is performed 
by increasing both [M/H] and [a/Fe] i n these models^ using 
Eq.[4| with A = 0.75 and A = 0.94 (see [Thomas et ^l200.3 ) 
for models and the TMJll models, respectively. Our mod¬ 
els, the TMJll models, and the CvD12 models all have the 
same native (MILES) spectral resolution, while the S07 line- 
strengths are given at the Lick/IDS resolution (see Table 7 
of S07F1- For a given age, metallicity, and [a/Fe], we have 
determined the behaviour of the S07 line-strengths at the 
MILES resolution, by smoothing our models to match the 
Lick/IDS resolution for each spectral feature, and then es¬ 
timating the corresponding variation of its line-strength at 
this resolution. 

It is evident in Fig.[2T]that the arrows corresponding to 
the various models do not have the same starting point in the 
index-index diagrams. These significant offsets among the 
predictions of the various models at scaled-solar abundance 
are, however, expected. These offsets reflect different model 
ingredients (e.g. the isochrones) and different values for the 
solar metallicity (see Section o among other aspects of the 
model construction. For instance, our (Mg, Fe, and Balmer 
lines) model predictions for an SSP with [M/H] = 0.0, age of 
13.5 Gyr, and [a/Fe]= 0.0, would be matched by the CvD12 
models with a 3.5 Gyr younger age, and slightly lower 
0.05 dex) total metallicity. Such differences strongly sug¬ 
gest that relative ages are more robu s t than absolute ages in 
the models (e.g., GlbsonetalJ 199^ VazdekisetaJj 1 2001 bl : 


ISchiavon et al.iT2002l . Mendel. Proctor fc Forb^l2007l ). Note 
also that such differences decreases with decreasing age. 
However, a detailed discussion on the origin of these off¬ 
sets between the models is beyond the scope of this paper. 
In the following, we focus on the relative responses of the 
various models to variations of the [a/Fe] abundance ratios 
(i.e., on the arrow directions, not their absolute values). 

Fig. l21l shows that the Mg6 index increases with [a/Fe], 
reflecting its strong sensitivity to [Mg/Fe]. The degree of in¬ 
crease is larger, by ~ 0.2 A, in our models and in those of 
TMJll, relative to the S07 and CvD12 models. All the iron 
lines show a mild dependence on [a/Fe] at fixed [Fe/H], as 
these lines are believed to be mostly sensitiv e to Fe abun¬ 
dance (see, e.g., S07, but also the discussion in iBarbuv et al.l 
I2OO3I : iGoelho et akiuOPsI about the effect of electron donors 
such as Mg on the continuum). While this does seem to be 
the case, Fig. 1221 shows that differences exist among model 
predictions. All models agree that Fe 5335 decreases with 


But note that the stellar library used in the S 07 models is of 
higher resolution than the Lick/IDS lljone 
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Figure 22. Same as Fig. 1211 but plotting the line-strengths of representative Balmer lines (top) and IMF-sensitive feat ures (bottom) as 
a function of the total metallicity indicator [MgFe]’. Note that the IMF-sensitive line-strengths are not provided by the lSchiavonI (120071 ! 
models. 


[a/Fe] (the amount of decrease being slightly stronger for 
our, relative to other, models). The behavior of Fe5270 is 
more model dependent, as this index remains constant with 
[a/Fe] for the CvD12 models, increases for TMJll, and de¬ 
creases in our models and in those of S07. The Fe4383 de¬ 
creases with [a/Fe] in all models, with a milder variation 
for our models and TMJll, relative to the CvD12 and S07 
models. Notice that, because of a possible dependence on 
Mg abundance (iKorn et al.ll200^ . the Fe4383 index might 
not be as clean an Fe indicator as as the Fe 5270 and Fe 5335 
indices (see S07). Indeed, Fig. [21] indicates that this index 
does not behave differently than the other Fe lines. Finally, 
Fe5015 is expected to have some sensitivity to the abun- 
dance of a-elemen ts, in particular to [Mg/Fe] and [Ti/Fe] 
llKorn et al.l l2005l ). However, no significant dependence of 
this index on [a/Fe] is seen in our models, CvD12, or the 
S07 models, while Fe5015 actually increases with [a/Fe] in 
the TMJll models. 

Fig. 1221 which includes the results for a number of age- 
and IMF-sensitive indices, shows several interesting features: 

H/3 line - TMJll, S07 and our models point to little de¬ 
pendence of H/3 on [a/Fe] at fixed [Fe/H]. In contrast, 
the CvD12 models predict that H/3 increases significantly 
with a-enhancement. Although the origin of this discrep¬ 
ancy is unclear, one should notice that H/3 decreases with 
velocity dispersion in ETGs (see Section 18.21) . and at 
fi xed velocity dispe r sion, as seen in the stacked spectra 
ofiii; Barbara et al.l ll2013l) . H/3 also mildly decreases with 
[Mg/Fe], by less than ~ 0.1 dex 0, from [Mg/Fe] ~ 0.1 to 
[Mg/Fe]~ 0.3. This seems to point to a genuine mild depen- 

We have estimated this amount of variation using the [Mg/Fe] 


dence of this index on enhancement, as predicted by S07, 
TMJll, and our models. We return to this issue in Sec¬ 
tion [8]2]3| 

Higher order Balmer lines - All the models agree in that 
HJf should increase with [a/Fe], although the amount of 
increase is larger for CvD12 and TMJll than for our models 
and S07. For H 7 F, TMJll predict that this index should 
have little dependence on [a/Fe], while our, S07, and the 
CvD12 models predict an increase of H 7 F with [a/Fe]. 

TiO IMF-sensitive features - We compare in Fig. l22l some 
representative IMF-sensitive features, i.e. NaD, TiOl, and 
Ti02, for which we can compare three different models 
(TMJll, CvD12, and ours). A compariso n between the 
CvD 1 2 models and the base-models from IVazdekis et al.l 
1I2OI2I) for ot her IMF-sensitive fea tures has recently been 
performed bv ISpiniello et aP (120141 ). The bottom panels in 
Fig. [22] show that all models agree on the prediction that 
TiO features should increase with [a/Fe]. This holds at ei¬ 
ther fixed [M/H] or fixed [Fe/H], as TiO bands are quite 
insensitive to total metallicity for old ages (see magenta ar¬ 
rows in the Figure, showing the effect of varying [M/H] by 
+0.3 at fixed age). Also, both the CvD12 and our mod¬ 
els agree that Ti02 should be less sensitive to [a/Fe], (when 
comparing the response to [a/Fe] with that to age, i.e. black 
and blue arrows in the Figure) than TiOl, making it a bet¬ 
ter IMF indicator than TiOl. The reader should be warned, 
however, that neither TiOl nor Ti02 show any significant 
variatio n with [Mg/Fe] in ETG s at fixed velocity disper¬ 
sion (see lLa Barbera et al ]B) . Interestingly, although not 
plotted here, we find that the sensitivity of Ti02 to [a/Fe] is 

stacked spectra of ETGs, for the two h ighest velocity dispersion 
bins, defined in lLa Barbera et al. I l2013l. see their Figure 7. 
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significantly smaller when a bottom-heavy IMF is adopted 
for computing the models. This is the case both for our mod¬ 
els and also the CvD12 models. 

NaD IMF-sensitive feature - This feature is found to de¬ 
crease with [a/Fe] at fixed [Fe/H] in agreement with the 
CvD12 models. In our models, NaD also decreases with 
[a/Fe] at fixed total metallicity, while at high [a/Fe] the 
increase of NaD with [M/H] is milder than at scaled-solar 
(see top-left panel of Fig. I17II . The net effect is that NaD 
decreases with [a/Fe], at either fixed [Fe/H] or [M/H]. It 
has been recently pointed out that an increase of total 
metallicity and [Na/Fe] (by -1-0.3 dex) might be able to 
explai n the high NaD l ine-strengths observed in massive 
ETGs Jjeong et al.ll20l3) . with no need for a bot t om-h eavier 
IMF in these galaxies. As noted bv ijeong et al.l (l2013l l. this 
result is based on the predictions from the TMJll models. 
Our new models support the prediction that NaD decreases 
with a-enhancement as in the CvD12 models. This suggests 
that it is necessary to invoke another source (e.g. IMF) in or¬ 
der to increase the model line-strengths sufficiently to match 
the high NaD of massive ETGs. However some caveats 
should be noted. For example, the NaD feature can be 
contaminated by interstellar absorption, although we note 
here that for low velocity dispersion ETGs with a Kroupa- 
like IMF (and selected to have low internal dust) our base 
model predictions for the Na D match exactly the obs erved 
strengths for this index (see (iLa Barbera et ani2013lH . In 
addition, we find that NaD increases, rather than decreases, 
with [Mg/Fe] at fixed velocity dispersion in massive ETGs. 
This is possibly due to a correlated va riation of [Na/Fe] 
and [ Mg/Fe] abundances in these systems dLa Barbera et ^ 
l2013h . 

In summary, Fig. l22l shows that although differences re¬ 
main among predictions from different models on the depen¬ 
dence of the Balmer lines on [a/Fe], our new models agree 
well with the expectations from most models compared here, 
i.e. that H/3 remains constant (as in S07 and TMJll), H 7 F 
increases (as in GvD12 and S07), and HJp increases (for 
all models) with [a/Fe]. Also, our models agree remarkably 
well with CvD12 models on the expected sensitivity of com¬ 
monly used IMF-sensitive features (TiOl, Ti02, and NaD) 
to [a/Fe]. 


8 COMPARISON WITH GLOBULAR 

CLUSTERS AND EARLY-TYPE GALAXIES 

8.1 Galactic Globular Clusters 

In Fig. [ 23 ] we compare selected line-strength indices mea¬ 
sured from our scaled-solar and a-enhanced model SEDs 
wi th those measur e d from the Galactic GC (GG G) libraries 
of ISchiavon et al.l (l2005l l and iKim et ^ (l2014h . Both the 
model SEDs and GGC data have been smoothed in order to 
achieve a common LIS resolution of 5 A (FWHM). 

The use of GGGs to test SSPs is particular useful; 
the majority of the GGCs in these two libraries have well- 
defined CMD-derived ages, and metallicities ([Fe/H]) and 
light element abundances such as [Mg/Fe], [G/Fe], [N/Fe] 
and [Ca/Fe] obtained from high-resolution spectroscopy (see 
e.g. iRoediger et al.l |2014| for a useful compilation of these 
quantities). 


For all the indices shown in Fig. 1231 at low metallic¬ 
ities ([M/H]< —1.0) the a-enhanced and base models are 
very similar. This is to be expected since the uncorrected 
MILES library stars for the base models in this regime are 
a-enhanced to a level similar to that of the a-enhanced mod¬ 
els. Additionally, at the lowest metallicities the index predic¬ 
tions for the scaled-solar models also converge towards the 
base- and a-enhanced models. Such behaviour is not surpris¬ 
ing as the sensitivity to the a-enhancement decreases with 
decreasing metallicity (see Fig. [Uj. This shows that the use 
of base models is also appropriate for studying metal-poor 
GGGs. 

The typical [M g/Fe] for GGC stars is ~ -1-0.4 

llRoediger et al.ll2014l l and this is in good agreement with 
location of the majority of the GGCs in the Mgb- (Fe) panel 
of Fig. 1231 Most of the GGCs lie on the model [Mg/Fe] = 
+0.4 locus. A couple of the most metal-rich GGCs (e.g. 
NGC6528 and NGC6553) lie between the scaled-solar and 
a-enhanced grids suggesting [Mg/Fe]+0.2 —1-0.3 which 
is consistent wi th the high-resolution [Mg/Fe] measurem ents 
of +0.25+0. llj|Caxret4a_etaJ j200ll : ICoelho et akll^OOlI l and 
+0.32+0.08 jBarbuv et al.lll992l ~) respectively. Our result for 
these tw o clusters are also in good agreement wit h i Lee et al.l 
(l2009blj (see their figure 13). 

[Ca/ Fe] is also enh a nced in the majority of GGCs 
(-^ +0.3; IRoediger et al.l (l2014| jj and this reflected in the 
location of GGCs in the Ca4227 - (Fe) panel of Fig. 1231 
[Ca/Fe] in GGCs is typically ~ 0.1 dex lower than [Mg/Fe] 
for a given cluster and this is well reproduced in the fig¬ 
ure. The fact that [Mg/Fe] and [Ca/Fe] in GGCs behave 
similarly when compared to the a-enhanced grids is both 
reassuring and not surprising - bo th Mg and Ca are a - 
elements. However, as pointed out bv iRoediger et akl (l2014li . 
the lack of any correlation between [Mg/Fe] and [Ca/Fe] 
in GGCs seems to suggest that, while both are indeed a- 
elements, they are possibly produced in different nucleosyn- 


track [Mg/Fel in massive galaxies fe.g..lVazdekis et al.l 

1997 

Trager et al. 

1 19981: ICenarro et al.ll2003l: iThomas et al.l 

2003 

Conrov 12013 

)• 


In terms of other indices, CN 2 (and CNi -not shown) 
is clearly stronger in the GGCs than the models, and the in- 


It is well established that GGCs (e 

g. Burstein et al.ll 19841: 

Cannon et al.l 19981: Roediger et al. 

l2014h and extragalac- 

tic GCs (e.g. iBurstein et al.lll984 

Brodie & Huchral 199ll: 

Kuntschner et al.l 2002a: Strader et al.l 2003: iBeaslev et al.l 

2 OO 4 I: ICenarro et al.l l2007bl: ISchiavon et al.l l2013l) show 


enhanced CN when compared to solar neighbourhood 
stars. Fig. JSl provid es suppo rt for the early re sults of 
I Lee fc Wort he vl (120051 ) (see also lBurstein et ^ri984h . which 
found a significant discrepancy when comparing Ca4227, 
CNi, CN 2 and NaD index strength from their a-enhanced 
SSP models and data for GCs in the Milky Way and in 
M31. This shortcoming of the current generation of a- 
enhanced SSP models can be understood when taking prop¬ 
erly into account the pheno menon of Multiple Populations 


(MPs) in Galactic GCs (see Gratton. Carretta fc Bragaglial 
l2012l|Piotto et al.l I2OI2I : I Cassisi fc Salarid I2OI3I . and refer- 
ences therein for detailed reviews on this topic). 

Essentially, the MP phenomenon is related to the pres¬ 
ence in each individual GC of two (almost coeval) stellar 
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Figure 23. Comparison of selected line-strength indices measured from our scaled - solar (green), a-en hanced (red ) and base (black) 
model SEDs with those measured from the Galactic GC libraries of ISchiavon et ^ (l2005l l (circles) and iKim et ahl ||2014| ~) (triangles). 
Both the model SEDs and GC data have been smoothed to match a common LIS resolution of 5 A (FWHM). Models with metallicities 
—2.27 sg [M/H] ^ - h0.06, ages 8 ^ t ig 14 Gyr and Kroupa Universal IM F are shown. Clusters are colour-cod ed by fMg/F e] taken from 
the compilation ofiRoediger et al ] ll2014h supplemented with data fror n l.lohnson et ahl |20^) ING C 62051. iBe hH |200,'jl (NGC6341), 
iKacharov. Koch fc Mcwilhami ll2013l l INGC 68641. [Rodgers &: Hardind lll99dl ('fCa/Fel: NGC 69811. [Kraft et ahl lIlOOsI l INGC 7006). 


generations, the first one born with a canonical a-enhanced 
heavy element distribution, the second generation charac¬ 
terised by anti-correlated CN-ONa abundance patterns, and 
by moderate (and - only - in a few cases very large) He 
enhancement. Due to the impact on the SED of molecu¬ 
lar bands associate d with both C and N (see the detailed 
discu ssion made by ISbordone et al.l l201ll and ICassisi et al.l 
l2013l) . one might anticipate that the light-element (anti- 
) correlations characteristic of the second generation stars 
in Galactic GCs will affect the strength of specific in- 
dic es such as GNi and GN 2 . This i ssue has been explored 
bv iGoelho. Percival fc SalarisI ll201lj) given certain assump- 
tion^T about the chemical abundance patterns of second 

They adopted extreme -although still realistic - light-element 
chemical patterns in order to maximise the possible expected im¬ 
pact of these chemical peculiarities on SSP modes. 


generation stars. These authors found that - regardless of 
the exact value of the He enhancement - as a consequence of 
the stronger GN and NH molecular bands in the spectra of 
second generation stars CN indices are significantly affected. 
Variations of the order of ~ 0.08 mag were seen between the 
CN indices predicted for a pure first generation stellar pop¬ 
ulation and those predicted for population of pure second 
generation stars (with extreme values for the light-element 
anti-correlations). 

Finally, unlike CN indices, the C24668 index of the 
GGCs is in good agreement with the model grids. As evi¬ 
denced by the superposition of the scaled-solar, a-enhanced 
and base models, it is clear that this index has little sensi¬ 
tivity to the a-abundance. 

In summary, we conclude that the LIS Mgb and iron 
indices of the GGCs are well-reproduced by the a-enhanced 
MILES models. Good agreement is also found for Ca 4227 in 
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GGCs. The model CN indices are largely insensitive to the 
a-abundance, and fail to reproduce the high GN values seen 
for GGGs at a given metallicity. The model C 2 4668 index 
is also insensitive to the a-abundance, but unlike GN repro¬ 
duces the behaviour of the GGC C 2 4668 measurements. 


8.2 Early-Type Galaxies 

8.2.1 Early-Type Galaxy Colours 

A strong constraint on the reliability of our models comes 
from the optical colours of nearby ETGs. Gurrent stellar 
population models, both based on theoretical stellar spectra 
and empirical libraries, still have difficulties in reproduc¬ 
ing the optical colour distribution of ETGs, with theoretical 
models showing the strongest discrepancies dMaraston et al.l 
I2OO9I : IConrov fc Gunnlboid : iRicciardelli et aLll2012ll . 

In Fig. [23] we present the colour distribution in the 
SPSS bands of a samp le of z~0.04 ETGs, as selected in 
iRicciardelli et ahl (I2OI2I '). These data are compared with 
our scaled-solar and a-enhanced SSP models for a range 
of metallicities and ages assuming a Kroupa Universal IMF. 
The scaled-solar models ( upper panels) hav e the s ame dif¬ 
ficulties as highlighted in iRicciardelli et aP (]2012f l for the 
base-models. This behaviour is expected since the scaled- 
solar models are very similar to the base models in the 
high metallicity regime required to fit luminous ETGs (see 
Fig-O, a region where the MILES stars have approximately 
scaled-solar abundance ratios (see Fig. [H]). The colours for 
the scaled-solar models are too red in g—r and too blue in 
u—g and r—i. This is a problem that cannot be solved by 
superimposing a young stellar population on top of the main 
population, since this would shift the whole set of colours 
to the blue and thereby worsen the mismatch with the u—g 
and r—i colours. 

As shown in Section lOl the a-enhanced models are 
bluer in the blue bands, and redder in the red, as seems to 
be required by the data. Indeed, the comparison with the a- 
enhanced models (lower panels of Fig. 1241) show a much bet¬ 
ter agreement with data in the g—r/r—i plane. However, the 
shift to the blue in u—g appears too dramatic to reconcile the 
model with these data, suggesting that an a-enhancement of 
[a/Fe]= 4-0.4 (as in the models presented here) is too high 
when compared to the observed colours of ETGs. However, 
the best solutions obtained with these models are more con¬ 
sistent than those obtained with the scaled-solar models. It 
is also worth recalling here that our models with varying 
abundance ratio do not extend to the blue end of the u fil¬ 
ter response, where we used our NGSL-based base model 
predictions. The a-enhanced models point to old ages in 
the colour-colour diagrams, whereas the scaled-solar models 
require significantly younger ages for the u—g versus g—r 
diagram. 

To explore the effect of more complex star formation 
histories, we perform a set of Monte Carlo (MC) simula¬ 
tions, where we allow all the stellar population parameters 
to vary: age, metallicity and IMF slope (using the bimodal 
parameterization). We adopt three-burst star formation his¬ 
tories, where the mass fraction involved in each burst is 
also a free parameter. The resulting predictions for 10000 
MC simulations are shown in Fig. [25] for scaled-solar mod¬ 
els (upper panels), for a-enhanced models (middle panels) 


and for a random combination of the two (lower panels). 
The mixture of scaled-solar and a-enhanced models gives 
better results as the solutions obtained for the two colour- 
colour diagrams are more consistent each other. The im¬ 
provement in the u—g/g—r plane with respect to the case 
where only SSPs are used, is mainly due to the contribu¬ 
tion of young stars and metal-poor populations. Neverthe¬ 
less, the models have difficulties in reproducing the part of 
the diagram lying at high u—g. The match in the r—i/g—r 
plane also improves significantly. Indeed, we see a clear red¬ 
dening of r—i, that is due to both the use of a-enhanced 
models and to varying the IMF slope (see Section 6.1 in 
iRicciardelli et aDi2012] l. The obtained solutions are consis¬ 
tent with recent studies based on absorption line-strengths, 
which find steeper IMF slopes [a/ Fe] values around ~ 0.2 
and old ages for massive ETGs (e.g..iLa Barbera et al ]i2013l : 
iSpiniello et al.ii2014] : iConrov &: van Dokkuin 2012bl L 


8.2.2 Full Spectrum Fitting 

The solar-scaled and a-enhanced model SEDs can be used 
in full spectral fitting techniques that can, in principle, im¬ 
prove the overall fits to data and be used to derive (in 
particular) the [Mg/Fe] overabundance observed in mas¬ 
sive galaxies. We have attempted to follow this approach 
by using the penalised pixel- fitting routine developed by 
iCappellari fc Emsellemi |]2004l '). We carried out this exer¬ 
cise o n data for Virgo galaxies published in lYamada et al.1 
(]2006tl . We perform our spectral fits on two ETGs, one rela¬ 
tively low-mass (NGC4489), and the other a giant elliptical 
(NGC4365) which have scaled-solar and [Mg/Fe]-enhanced 
abundances, respectively. 

At first we tried to perform spectral fits over large wave¬ 
length ranges (e.g. more than 1000 A), but we were unable to 
reproduce the [Mg/Fe] values a nd trends within galaxies re¬ 
ported bv lYamada et all l]2006li . This is not surprising given 
the fact that the information contained in relevant Mg- and 
Fe-sensitiv e features gets dilute d with increasing wavelength 
range (e.g. IWalcher et al.11200^ . We therefore turned to fit 
only those regions where we expected the effect to be most 
noticeable (i.e. the Mg6 - Fe5335 range; see for example 
Fig. [ini and Fig. [TT]) . We therefore performed the fit within 
the range AA 4800 — 5500 A. The results of this exercise are 
shown in Figs. 1261 and l27l Each figure shows our attempts to 
fit the observed spectrum using (i) scaled-solar models only, 
(ii) a-enhanced models only, and (iii) a combination of the 
two using the same age and metallicity for the two popula¬ 
tions. In Fig. [2Hwe see that, for NGC4489 - a low velocity 
dispersion galaxy - pure scaled-solar and the mixed mod¬ 
els reproduce the spectrum fairly well, although pPXF has 
some preference for the true scaled-solar models. Fits with 
a-enhanced models only cannot match all the observed ab¬ 
sorption features. The case of the massive galaxy NGC4365 
{a — 245 kms“^) is different (see Fig. ilZ]). Here a-enhanced 
models are preferred over the scaled-solar models. It is im¬ 
portant to note that the higher velocity dispersion of this 
galaxy makes it harder to identify any mismatches with the 
templates, as can been seen in the residual spectra. 

The overall conclusion from this exercise is that our 
models together with full-spectral fitting can be used to es¬ 
timate the [Mg/Fe] enhancement. However the spectral fits 
are sensitive to the wavelength range employed to estimate 
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Figure 24. Observed colour distribution in the SDSS bands of a sample of ETGs at 0.04. The contours enclose 20, 40, 60 and 95 
per cent of the galaxies. The cross indicates the location of the median colours. The synthetic colours derived from scaled-solar models 
(upper panels) and Q-enhanced models for a Kroupa Universal IMF are overplotted. Only SSP models with ages in the range 3—13.5 Gyr 
and metallicities: —0.35, -1-0.06, -hO.26 have been considered. 


the a-enhancement. Short wavelength regions around the 
main spectral features (e.g. Mgb, Fe 527 0) are preferred over 
long baselines (see also the discussion in IWalcher et al.ll2009l 
and ICezari 3 120131 on the dependence of spectral fitting re¬ 
sults with the wavelength range). 

8.2.3 Line indices 

Fig. [28] shows three classical index-index diagrams with 
the predictions of the three models, base, scaled-solar and 


a-enhanced. Overplotted are the index measurements in 
the LIS- 8 . 4 A system for the e arly-type galaxy sample of 
ISanchez-Blazauez et al.l (I2006al 1. The model grids obtained 
when plotting H/3o against Mgb, (Fe) and [MgFe] indices, 
which are sensitive to Mg, Fe and total metallicity, respec¬ 
tively, provide similar age/metallicity solutions for the most 
massive galaxies, i.e. those galaxies with the largest abun¬ 
dance ratios and largest [MgFe] index values. 

We see that most of the massive galaxies fall below the 
horizontal line corresponding to the oldest stellar popula- 
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Figure 25. The same ETG colour distributions of Fig. I24l are compared with the predictions from composite stellar population models. 
The black contours enclose 20, 40, 60 and 95 per cent of the MC simulations of multiple bursts models, where age, metallicity, IMF 
(bimodal) slope and strength of the bursts are allowed to vary. The upper panels show the predictions for the scaled-solar model, the 
middle panels display those of o-enhanced models, whereas the bottom panels show the predictions from MC simulations where a random 
combination of the two set of models is considered. 
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Figure 26. Spectral fitting of the low velocity-dispersion galaxy 
NGC4489. The upper panel shows the fit obtained with scaled- 
solar models. The middle panel shows the fit obtained with the a- 
enhanced models with [Mg/Fe]= -1-0.4. The bottom panel shows 
the fit with both set of models included in the database. The 
residuals are shown (in green), keeping the same scale. 

tion (14Gyr) of the scaled-solar model grids. The same ap¬ 
plies to the base model grids, which at high metallicities 
have an approximately scaled-solar [Mg/Fe] abundance ra¬ 
tio. However, the locus of these galaxies falls in between the 
lines corresponding to the oldest SSPs of the scaled-solar 
and the a-enhanced models, consistent with old ages, to¬ 
tal metallicity above solar and an abundance ratio around 
[Mg/Fe]~ -1-0.2. Note that we only employed models with a 
Kroupa Universal IMF for the grids of Fig.[2H]as our present 
purpose is not fit these galaxies, but rather to illustrate the 
advantages of the new models presented here. We could have 
obtained slightly younger ages for the most massive ETGs 
had we employed grids with steeper IMF slopes, according 
to the relation between the slope and the velocity dispersion 
according to (Ferreras et al.ll20F^ : iLa Barbera et al.l 20131 : 
ISpiniello et al. 2014h (see also the results we show and dis¬ 
cuss shortly in Fig. 1^ . 

It also is worth noting in Fig. [28] (top panel) that the 
scaled-solar line that corresponds to solar metallicity, i.e. the 
third vertical line starting from the right, is very similar to 


Figure 27. Spectral fitting of the giant elliptical galaxy 
NGC4365. See Fig. 1261 for the details on the panels. 


that corresponding to the base model. Similar behaviour is 
seen for the more metal-rich lines with [M/H]= -1-0.15 and 
-1-0.26, although in these cases the lines corresponding to 
the base models are slightly shifted toward the solar metal¬ 
licity value with respect to the scaled-solar models. We at¬ 
tribute this small difference to the fact that many of the 
MILES stars have [Mg/Fe] v alues slightly below scaled-solar 
at these high metallicities (|Milone et al.l[201ll) . For lower 
metallicities, the vertical lines corresponding to the base 
models are shifted toward a higher metallicity with respect 
to the two, the scaled-solar and the a-enhanced model lines 
(top panel. Fig. 1281) . Notice that for the base models the 
stars are selected according to their [Fe/H[ values. However 
at these metallicities the MILES stars are over-enhanced in 
the [Mg/Fe[ abundance. Therefore the stars have a total 
metallicity that is higher than that of the stars employed in 
the scaled-solar and a-enhanced models according to Eq. (4] 
Fig. [29] is similar to Fig. [28] but shows representative 
Balmer lines, H/3o, H/?, HJp, and H 7 F (panels a-d), as a func¬ 
tion of the total metallicity indicator [MgFe]'. The grids with 
dotted linetype in the Figure show models with a bottom- 
heavy, rather than Kroupa-like, IMF. Green and red grids 
correspond to scaled-solar and a-enhanced models, respec- 
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Figure 28. Various Mg and Fe index-index diagrams for 
three different model sets; the base models (black), the 
scaled-solar models (green) and the a-enhanced models (red), 
with [Mg/Fe] = -1-0.4. Overplotted are the line-strength 

indices measured in the sample of ISanchez-Blazauez et al] 
i( 20063)- The metallicity increases from left to right [M/H] = 
—0.35, —0.25, -h0.06, -1-0.15, -1-0.26, whereas the age increases from 
top to bottom 2.0, 2.75, 5.0, 7.5,10.0,12.5,14.0 Gyr. All the mod¬ 
els have a Kroupa Universal IMF and have been smoothed to 
match the LIS-8.4A. 


lively. In each panel, we overplot data-poin t s for 18 stacked 
spectra of ETGs from iLa Barbara et al] (|2013l) . ranging 
from a galaxy velocity dispersion of 100 to ~ 300 kms“^. 
Both model and data line-strengths are plotted at the same 
spectral resolution of 200 kms“^, added in quadrature to 
the native resolution of the SDSS spectrograph (~ 2.8 A at 
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Figure 29. Line-strengths of Balmer lines, H/3, H/3o, and 

Hyp are plotted as a function of the total metallicity indica¬ 
tor [MgFe]', for scaled-solar (green) and o-enhanced SSP mod¬ 
els (red). Light-solid and dark-dotted grids correspond to mod¬ 
els with standard (bimodal; F^ = 1.3) and bottom-heavy (bi- 
modal; F;, = 2.8) IMF, respectively. The grids show the effect 
of changing total metallicity (moving from left to right in each 
panel) and age (increasing from top to bottom in each panel), 
respectively, as labelled. Blue dots, connected by a blue line, 
mar k the line-strengths for 18 SDSS stacked spectra of ETGs, 
from iLa Barbera et al.l l|2Q13h . in the range of galaxy velocity dis¬ 
persion from ~ 100 to ~ 300 kms“^. Line-strengths have been 
corrected for all models and spectra to the same <j of 200 kms~ ^ 
(i.e. the same resolution employed in iLa Barbera et ^ 1 120131) . 
summed up in quadrature to the SDSS spectral resolution. 


A ~ 5000 A). The comparison of models and data can be 
summarised as follows. 

i. The model H/3 and H/3o indices behave differently with 
a-enhancement, in that the former is nearly independent 
of [a/Fe], while H/3o decreases with a-enhancement. As a 
result, one can explain the observed H/3o of our stacked 
spectra of ETGs with a-enhanced models, while for H/3 
a bottom-heavy IMF (see dotted grids in panel b of the 
Figure) seems to be required to account for the observed 
line-strengths. However, as noted above (Sec. [Tj, the be¬ 
haviour of H/3 with [a/Fe] is also quite model dependent, 
as the CvD12 models predict that this index will increase, 
ra ther than remain const a nt, wi th [a/Fe]. Also, as reported 
in ICervantes &: Vazdekij ll2009l) . different theoretical stel¬ 
lar libraries lead to different predictio ns for the response o f 
H/3 and H/3o to [a/Fe]. Irideed, u sing iMunari et alJ ll2005l) . 
rather than ICoelho et ahl (l2005l ) stellar models, one would 
predict the opposite [a/Fe] dependen ce of H/3 and H/3o than 
that seen in Fig. (see section 4 of ICervantes fc VazdekisI 
I2OO9I) . These differences originate, at least in part, in the 
residual seen toward the red end of the feature bandpass 
of the H/3 and H/3o index definitions, as illustrated in the 
flux-ratio of Fig. 1161 This residual corresponds to a metallic 
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contamination within the red wing of the H/3 feature, which 
includes Fe-peak elements such as chromium (see also Sec¬ 
tion lOl) . The discrepancies betw een the predict i ons fo r the 
H/? and H/3o lines based on the iMunari et alJ ll2005l ) and 
ICoelho et alJ ll2005l l stellar spectra likely stem from differ¬ 
ences in the adopted line lists employed for the computation 
of these two sets of model spectra, and are presently being 
investigated (Sansom & Coelho, in prep.)- 

ii. The position of the stacked spectra with respect to 
scaled-solar model grids is approximately the same in both 
the H/3 vs. [MgFe]'and H/3o vs. [MgFe]' diagrams. This, to¬ 
gether with the fact that both H/3 and H/3o exhibit a simi¬ 
lar variation with [a/Fe] at fixed velocity dispersion in the 
stacked spectra of ETGs (see Sec. 0, suggests that both 
H/3 and H/3o have more similar dependencies on [a/Fe] than 
suggested by the models. 

hi. The H5 f index behaves consistently with H/3o in the 
sense that both indices, in the stacked spectra of ETGs, 
could be explained by a metal-rich, old and significantly q- 
enhanced ([a/Fe] ~ -1-0.3 dex) SSP. 

iv. In contrast, the behaviour of H 7 F with [a/Fe] seems 
to be too extreme in the models. The observed H 7 F of 
the stacked spectra would require [a/Fe]~ 0.2 dex for an 
old (~ 14Gyr) and metal-rich ([M/H] ~ -1-0.3) SSP, in 
contrast to the higher [a/Fe] implied by H/3o and H5 f. 
The position of the stacked spectra with respect to the 
scaled-solar H 7 F grid suggests that the “true” dependence 
of H 7 F on [a/Fe] should indeed be mild. Note however 
that the blue pseudo-continuum of the H 7 F line is located 
on a strong CH absorption, i.e. the G-band at ~ 4300 A. 
Therefore, non-solar carbon (and nitro gen) abundance ra¬ 
tios reported in massive galaxies (e.g. iGonrov et al.l I20l4 
IWorthev. Haitia n fc Servenll20 1^. which might depend on 
environmen t I Sanchez-Blazauez et al.ll2003l : ICarretero et al.l 
l2004l200'itl . may have a noticeable effect on measured H 7 F 
index strength. 


In summary, Fig.l^Hlpoints to the well-known issue that 
SSP model predictions for different Balmer lines (and in¬ 
deed also different definitions of the same line, e.g. H/3 and 
H/3o) provide different estimat es of the age of t he old stellar 
populations of ETGs (see also ISchiavon et al.l l^2004^ for an 
alternative view, who argues that the disagreement between 
high order Balmer lines might be related to the presence 
of young stellar populations.) The Figure shows that nei¬ 
ther the a-enhanced models nor models with a varying IMF 
seem able to solve this issue. Some clues come from the com¬ 
parison of the [a/Fe] dependence of Balmer lines among dif¬ 
ferent models (see Fig. I22II and the observed line-strengths 
of ETGs, as reported above. However, these clues should be 
only interpreted in a qualitative manner and a more quan¬ 
titative comparison of models and data is desirable. In fact, 
the discrepancies between data and models in Fig. 1291 might 
occur due to different effects, such as a variation of abun¬ 
dances ratios of specific chemical species in the stellar atmo¬ 
spheres, and/or more complex star-formation histories than 
pure SSPs. Also, the observed differences in the predictions 
from the various models suggest differences in the theoretical 
stellar spectra that are employed to obtain the differential 
corrections, which might come in part from the adopted line 
lists. 


8 . 2.4 The proxy for [Mg/Fe] 

As discussed above (Section 18. 2. 211 . our new a-enhanced 
models can be used to estimate the [Mg/Fe] of ETGs 
through spectral fitting, provided one defines carefully the 
spectral range to fit. In this section, we discuss alternative 
methods to estimate [Mg/Fe], by either comparing the line- 
strengths of selected (Mg and Fe) spectral features to pre¬ 
dictions of e ither a-enhanced (as in , e.g., TMJll) or base 
models fsee iLa Barbera et al.l 120131 /. To this end, we rely 
on a subsample of 4,671 nearby {z ~ 0.07) ETGs from 
the SDSS-based SPIDER sample, restricting the analysis 
to the best-quality spectra, with S/N > 25 per A in the 
spe ctral regions of the Mg b, Fe5270 and Fe5335 features 
fsee lLa Barbera et al ][ml and references therein). For each 
spectrum, we obtain different estimates of [Mg/Fe] as fol¬ 
lows: 

“True” [a/Fe] - We fit [MgFe]', Mgfo and Fej^ line- 
strengths with either TMJll or our a-enhanced models, fix¬ 
ing the age of the SSP models in the ht to the luminosity- 
weighted age me asured with the softwar e STARLIGHT for 
each galaxy fsee lLa Barbera et al. I I 2 OI 3 I for details). 

Proxy for [a/Fe] - At fixed (STARLIGHT) age, we es¬ 
timate the metallicity by fitting alternately the Mg6 
and Fe3 line-strengths using either the base or the 
scaled-solar models, without the differential correction 
for the alpha-enhancement effect (i.e. [M/H]Mgb and 
[M/H]Fe 3 ). The difference of the two metallicity esti¬ 
mates, [.Z’Mo/.^Fe] = [M/H]Mgb —[M/H]Fe 3 . E the proxy for 
[Mg/FerT The computation is repeated for three sets of 
reference models, i.e. our scaled-solar, base models and 
base models constructed with the PadovaOO isochrones. 
We refer to these three proxies as [ZMg/ZFe]ss(BaSTi) > 
[ZMg/ZFe]Base(BaSTI) and [ZlVIg/ZFe]Base(PadovaOO) 7 respec¬ 
tively. 

As discussed in iLa Barbera et akl ll2013lf , since for an 
a-enhanced population one has [M/H]m 96 > [M/H]Fe 3 , es¬ 
timating the proxy may require some extrapolation of the 
reference models above the highest available [M/H] in the 
models. This may hamper the estimate of the proxy for ex¬ 
tremely high_[a/Fe]^_On the o ther hand, as shown in hg- 
ure 6 of iLa Barbera et al. liioTi) , the [ZMg/ZFe]Base(Padova00) 
proxy correlates remarkably well with the [a/Fe] estimated 
from models (TMJll) that take [a/Fe] explicitly into ac¬ 
count. This suggests that for several applications one does 
not need to use a-enhanced models to estimate the [Mg/Fe], 
provided that a suitable calibration of the proxy into the 
“true” [a/Fe] is perfo rmed. To this ef fe ct. F ig. 1501 ex¬ 
pands on figure 6 of iLa Barbera et ah] 1 I 2 OI 3 II , compar¬ 
ing the [a/Fe] estimate from our new models with the 
three proxies [ZMg/ZFe]sS(BaSTI), [ZMg/Zpe]Base(BaSTI)and 
[ZMg/ZFe]Base(PadovaOO) (panels a-c, respectively), and com¬ 
paring it to the [a/Fe] from the TMJll models 
([a/Fe]TMjii; panel d). Solid lines in the Figure show best- 
htting lines to the data, while dashed lines mark the ±1- 
sigma scatter around median trends. The Figure confirms 


1"* Fe3 = i(Fe4383 + Fe5270 -I- Fe5335); iKuntschneJ ||2000|~) . 

Note that we keep for the proxy th e same notation we 
employed in our pre vious works (e.g., lYamada et al.l l2006l: 
iLa Barbera et ^l2013l) . 
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Figure 30. The “true” [o/Fe], as estimated by our new models (see the text), is compared to different proxies for [o/Fe] (panels a—c), 
and the [a/Fe] estimated with the TMJll models (panel d). We compare [a/Fe] to the proxy for scaled-solar (a), base-BaSTI (b) and 
base-Padova (c) models. For each panel, grey dots correspond to a sample of 4, 671 ETGs selected from the SPIDER sample (see the 
text), solid curves are the best-fitting lines to the data (using [a/Fe] as dependent variable in the fitting), dashed curves are the ztl-cr 
scatter around the median trends, while horizontal and vertical dashed lines mark the values of zero on the x- and y-axes. The scatter 
in [a/Fe] around the best-fitting lines amounts to 0.025, 0.05, 0.04 and 0.025 dex for panels a to d, respectively. 


the remarkable correlation of the true [a/Fe] with the prox¬ 
ies. However, differences do exist among different proxy es¬ 
timates, in that (i) the slope and zero-point of the best¬ 
fitting line are slightly different from one proxy to another, 
the former ranging from ~ 0.56 to ~ 0.62^H and the latter 
from ~ 0 to ~ 0.05. Notice also that the scatter around the 
proxy - [a/Fe] best-fitting lines is smallest for our scaled- 
solar models (-^ 0.025 dex in [a/Fe]), and increases to 
~ 0.04 dex for [ZMg/ZFe]Base(Padova 00 ), and to ~ 0.05 dex 
for [ZMg/ZFe]Base(BaSTi)- The Scatter is largest for the base- 
BaSTI models because of a prominent tail of objects having 
too low [ZMg/ZFe]Base(BaSTi) at given [a/Fe] in panel b (rel¬ 
ative to panel a) of Fig. 1301 The origin of this tail can be 
understood through Fig. 1311 comparing the dependence of 
Mg 6 and Fe3 line-strengths on total metallicity for different 
non a-enhanced models. While the slope of the Fe3-[M/H] 
relation turns out to be approximately the same for all mod¬ 
els, the behaviour of Mg& vs. [M/H] is more model depen¬ 
dent, with the shallowest slope for the base-BaSTI models. 
As seen by comparing the green and black curves in the top 
panel of Fig. 1311 at high [M/H] (> 0) one obtains similar 
[M/H]m 96 and [M/H]Fe 3 estimates when using either base- 
BaSTI or scaled-solar models. In contrast, at low [M/H], 
the [M/H]m 96 estimate is always lower (for given Mg 6 line- 
strength) for base-BaSTI than scaled-solar models, leading 
to lower [ZMg/ZFeXs, explaining the tail of low proxy values 
in panel b, relative to panel a, of Fig. 1301 

Fig. 1301 panel d, also shows that a tight correlation ex¬ 
ists between the “true” [a/Fe] estimates from our models 
and from the TMJll a-enhanced models, with a scatter of 
only 0.025 dex in [a/Fe]TMjii. Notice that the slope of the 
best-fitting relation between the two estimates of [a/Fe] is 
larger than one (~ 1 - 2 ), implying that using our models 

Note that in figure 6 of iLa Barbara et al.l ll2013l l we reported 
a gradient of ~ 0.55, which is different from the values in Fig. 1301 
In that paper we compared the proxy to the TMJll [a/Fe] esti¬ 
mates, and did not apply the same S/N cut to select ETGs, as in 
the present work. 


leads to values of [a/Fe] larger than those from the TMJll 
models. 

In summary, we find that relative estimates of [Mg/Fe], 
with accuracy < 0.05 dex, may not require (for most applica¬ 
tions) stellar population models with enhanced [a/Fe], while 
absolute estimates of [a/Fe] require “truly” a-enhanced 
models. However, absolute estimates of [a/Fe] are affected 
by differences in the adopted stell ar populations mod els 
(consistent with what found, e.g., bv IConrov et al.ll2014^ ■ 


9 THE MILES WEBPAGE 

In an effort to facilitate the use of our models, we have in¬ 
tegrated these new predictions into our dedicated website 
(http://miles.iac.es). The main new feature is the ability 
to compute spectra for more complex SFHs, with varying 
IMF and a-enhancement. As for the current set of models, 
our website will deliver spectra at the desired instrumental 
setup decided by the user, as well as provide line-strength 
and colour predictions among other observables. 

9.1 Complex models with multiple bursts 

The complex SFHs used for Fig. [5S] have been obtained 
through the SFH facility implemented in the MILES web¬ 
site. This tool allows the user to retrieve the model spectrum 
corresponding to a given SFH, having an arbitrary num¬ 
ber of bursts, with differ ent ages, metallicities and burst- 
strengths (see Section 3 in iRicciardelli et al.ll20ld ). The SFH 
required by the user can either follow some parametric de¬ 
scription or be user-defined. For the parametric SFHs four 
different cases are allowed: multiple bursts (analogous to the 
user-defined case but limited to 5 bursts only), a truncated 
SFH, in which the star formation rate (SFR) is assumed 
to be constant between the time of formation and trunca¬ 
tion, an exponential SFH and an exponential plus bursts, 
where up to five bursts can be added on top of the expo¬ 
nential SFR. The user-defined SFHs are designed to handle 
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Figure 31. Behaviour of Mg6 (top) and Fe3 (bottom) as a func¬ 
tion of total metallicity for the base-BaSTI, scaled-solar, and 
base-Padova models, as labelled in the upper panel. Solid and 
dotted curves correspond to ages of 12 and 6 Gyr, respectively, 
typical for the spectra of ETGs used for Fig. I3UI 


complex SFHs, with an arbitrary number of bursts permit¬ 
ted. In the new web implementation, the user-defined SFH 
has been upgraded to include also variable IMF slope and 
a-enhancement. 

Given the mass fraction, the metallicity, the IMF slope 
and the a-enhancement at different ages, the tool computes 
the corresponding SED, Sx, according to: 

Sa = ^/^,SA(ti,[M/H].,[a/Fe]„-I>(ri),/.J (7) 

i=l 

where Nb is the number of input bursts, fm^ is the fraction 
of mass formed in the i-th burst, at age ti, with metallicity 
[M/H]^, IMF slope Fi (or Ft^ in the case of a bimodal IMF) 
and enhancement [a/FeJ^, which in the present study can be 
+0.0 and +0.4. SxlU, [M/H]„ [a/Fe]„-I>(ri),is the i-th 
SSP spectrum, which is computed with an isochrone, la, 
with similar a-enhancement. 


10 SUMMARY 

This paper is the second, and last, of a series 
where we present evolutionary stellar population syn- 
thesis models based on the empirical library MILES 
lISanchez-Blazauez et al.ll2006l : ICenarro et al.ll20073 ). which 


represents a significant improvement over previous em¬ 
pirical libraries. In the first paper (Paper I) we devel¬ 
oped ’’base models”, employin g the PadovaOO scaled-solar 
isochrones jGirardi et ahll^OOl) and MILES to predict spec¬ 
tra of intermediate- and old-aged stellar populations in the 
range AA 3540.5-7409.6 A, at moderately high resolution 
(FWHM=2.51 A). Following the abundance pattern of the 
solar neighborhood the computed model SEDs are enhanced 
in a-elements at low metallicity, whereas around solar metal¬ 
licity these predictions are approximately scaled-solar. The 
base models are therefore inconsistent at low metallicities. 
In the present work we further developed our models based 
on the MILES library, as well as theoretical stellar spectral 
libraries used in a differential way, to provide self-consistent 
predictions at all metallicities. We emphasise that our mod¬ 
els can be regarded as self-consistent in the sense that both 
ingredients, the theoretical stellar spectra and the isochrones 
are computed with the same overall a-enhancement. We 
summarise below these model developments and results: 


i) We compute d fully-theore ti cal m odels based on the 
stellar library of [Coelho_et_ay ([200^, with its extension 
to cool stars of Coelho et al.l ( 20071). to obtain SSPs- 
based differential spe ctral corrections ( Cervantes et al ]|2007l : 
IWalcher et al.l 20091). W e employed the BaSTI isochrones 
dPietrinferni et ^ l2004l . l2006l ). calculated with similar 
[a/Fe] abundance ratios to the stellar atmospheres, con¬ 
verted to the observational plane using exte nsive empiri¬ 
cal photom etric stellar libraries, mostly from lAlonso et al.l 
lll996l . [l9^ ). To improve the metallicity coverage around the 
solar value we computed a new set of BaSTI isochrones with 
metallicity Z = 0.024 ([M/H] = +0.15). The obtained theo¬ 
retical differential corrections are applied to the model SEDs 
computed with MILES, taking into account the [ Mg/Fe ] 
determinations of the library stars of iMilone et al.l ll201ll) . 
Therefore, at solar metallicity the scaled-solar models are 
truly scaled-solar and fully based on MILES, whereas for the 
a-enhanced models, we apply the differential corrections to 
the model SEDs computed with MILES. At low metallici¬ 
ties, the opposite is true; the MILES library is intrinsically 
a-enhanced and we correct the model SEDs using differen¬ 
tial corrections to obtain scaled-solar SEDs at low metal¬ 
licity. We therefore compute self-consistent models for two 
abundance ratio values: [a/Fe]= 0.0 and [a/Fe]= +0.4. 

ii) In addition to the scaled-solar and a-enhanced mod¬ 
els, and to the models we published in Paper I using the 
PadovaOO isochrones, we also computed base models em¬ 
ploying the BaSTI scaled-solar isochrones. 

iii) We have assessed the quality of the new models, estab¬ 
lishing the parameter space where the models are safe. The 
new models cover the metallicity range —2.3 ^ [M/H] ^ 
+0.4, although the model predictions, and in particular 
those of the scaled-solar and alpha-enhanced models, are 
generally safe within the range —1.49 ^ [M/H] ^ +0.26 
(within the constraints listed in Tableland shown in Fig. [6] 
and Fig. 0. Also, the spectral range blue-ward of 3800 A is 
unsafe for the scaled-solar models with sub-solar metallic¬ 
ities and for the a-enhanced models with solar and higher 
metallicities. We strongly suggest the interested user of these 
predictions to read the notes on the unsafe model parameter 
ranges that we provide in Table [T] and Table (2] It is particu- 
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larly important to do so when the obtained fits to the data 
fall within these unsafe ranges. 

iv) We provide SSP spectra for a suite of IMF shapes: 
single-power law (unimodal ), low-mass (< 0.5Mp,) tap ered 


single-power law (unimodal ), low-mass (< 0.5Mp,) tap 
IMF (bimodal) as defined in lVazdekis et all ( 199^ 200^ 


and 


the Kroupa Universal and Revised IMFs of Kroupal l ioOll) . 
For the unimodal and bimodal IMFs we have varied the 
logarithmic slope from 0.3 (top-heavy) to 3.3 (very bottom- 
heavy). A unimodal shape with slope 1.3 matches approx- 
imataley the Salpeter IMF, whereas a bimodal shape with 
the same slope is very close to the Kroupa Universal IMF. 
The predictions for the models adopting a unimodal IMF 
are safe for slopes smaller than 2.3, taking into account the 
above constraints. 


v) We have characterised the behaviour of the newly syn¬ 
thesised SSP model spectra. We find that, at constant total 
metallicity, [M/H], the a-enhanced models, which have lower 
iron content and therefore have lower opacities, show an ex¬ 
cess of flux blue-ward of ~ 4500 A with respect to the scaled- 
solar models. Besides this excess in flux, the a-enhanced and 
the scaled-solar models differ signihcantly in prominent fea¬ 
tures that are mainly related to the Call H-K and several 
CN and CH molecular absorptions. Both the blue flux ex¬ 
cess and the differences seen in these features, which mostly 
originate in the stellar atmospheres, increase with increasing 
age and increasing metallicity. In comparison to the effects 
of the stellar spectra the effects of the isochrones are of sec¬ 
ond order, and mostly impact on the continuum over the 
whole spe ctral range co v ered by the models, confirming the 
results by ICoelho et al.l (l2007l ) . 

vi) We confirm previous results that the combined [MgFe] 
and [MgFe]' are virtually insensitive to the [a/Fe] abun¬ 
dance ratio, and therefore are suitable for measuring the 
total metallicity. 

vii) We also confirm previous findings that the Hy and 
H(5 Balmer line indices are more sensitive to [a/Fe] than 
H/3. However, we show that such sensitivity is due to the 
pseudo-continua employed for the higher order Balmer index 
definitions, which are placed on relevant features within the 
blue flux excess of the a-enhanced models. The narrower the 
Balmer index definition, the lesser the sensitivity to [a/Fe]. 

viii) As expected, the Fe-dominated indices, such as 
Fe4383, Fe5015, Fe5270 and Fe5335 are weaker, and Mg6 
stronger in the a-enhanced models. This behaviour becomes 
more pronounced with increasing age and increasing metal¬ 
licity. Such trends do not apply to Mgi and Mg 2 , which are 
stronger in the scaled-solar models, as these indices show 
a greater sensitivity to carbon than [a/Fe]. These and all 
the other non-a elements are depleted with respect to the 
scaled-solar models. A similar behaviour is seen for C 2 4668 
and CN 2 (and CNi), which depend on C and N abundances 
respectively. Also in our models the NaD index behaves as 
the Fe-dominated indices. 


ix) We also see that at ages less than ~ 0.1 Gyr there is 
a signihcant increase in the strengths of various lines, most 
noticeably for the IMF-sensitive indices. This is as a conse¬ 
quence of the influence of supergiant stars at these young 
ages. Among the IMF indicators we find that, unlike TiOi, 
the Ti02 index is almost insensitive to [a/Fe], in good agree¬ 
ment with other authors. We also find that the Ti02 index, 
and other indices that show little sensitivity to age when 


employing a standard IMF slope (i.e 1.3), become more sen¬ 
sitive to this parameter for bottom-heavier IMFs. 

x) The base models, which approximately follow the abun¬ 
dance pattern of the Milky Way (MILES) stars, show in¬ 
dex values that are more similar to those of the a-enhanced 
models at low metallicities ([M/H] < —0.4), despite the fact 
that the employed isochrones are scaled-solar. On the other 
hand, at high metallicities the base models provide index val¬ 
ues that are very close to the scaled-solar ones. In general, 
the indices of the base models computed with BaSTI and 
PadovaOO isochrones are very similar, with the line-strengths 
obtained with BaSTI being slightly larger for the Balmer 
lines and slightly smaller for the metallicity indices. There¬ 
fore, for a given dataset, somewhat older ages and higher 
metallicities would be inferred from the BaSTI-based mod¬ 
els than from PadovaOO-based models. 

xi) A similar behaviour is seen for the colours, with the 
models based on BaSTI being bluer than those based on 
PadovaOO. Overall the agreement is good for the old stellar 
populations. For ages where the AGB contribution is more 
significant we find the larger discrepancies as a result of the 
varying descriptions for this evolutionary phase. 

xii) We use the base models as a reference for assessing the 
behaviour of the scaled-solar and a-enhanced colours. We 
hnd that the a-enhanced models are bluer in the blue, by 

0.2 mag in u—g and ~ 0.05 mag in g—r for old populations, 
and redder in the red, by ~ 0.02 mag in r—i, with respect to 
the base models. The best agreement between the base and 
a-enhanced models is reached at low metallicities. On the 
other hand, for the scaled-solar models the best agreement 
is reached at high metallicities. 

xiii) We compared our newly-synthesised SSP spectra and 
line-strengths with other models in the literature. In partic¬ 
ular we find that there is an overall agreement with the SSP 
spectra of the CvD12 models. We find that the overall ef¬ 
fect of increasing ]a/Fe] is similar to our models, although 
there are signihcant differences in the relative effect on dif¬ 
ferent spectral features. We also compared our line-strength 
indices with the predictions of the TMJll and S07 models. 
As expected we hnd signihcant offsets among the predic¬ 
tions of the various models, even at scaled-solar abundance. 
At hxed [Fe/H] the amount of increase of the Mg& index 
is larger, by ~ 0.2 A, for our models and the TMJll mod¬ 
els, relative to the S07 and CvD12 models. The iron lines 
show a mild dependence on [a/Fe], but we hnd that dif¬ 
ferences exist among the various models. Whereas all the 
models agree on Fe5335 and Fe4383, which decrease with 
[a/Fe] (although with varying sensitivities), the Fe5270 and 
Fe5015 indices are more model dependent. Although diher- 
ences remain among predictions from different models on 
the dependence of Balmer lines on [a/Fe], our new models 
agree well with expectations from most models compared 
here, i.e. that H/3 remains constant (as in S07 and TMJll), 
Hqp increases (as in GvD12 and S07), and HJp increases (for 
all models) with [a/Fe]. Finally our models agree remark¬ 
ably well with CvD12 models on the expected sensitivity 
of commonly used IMF-sensitive features (TiOl, Ti02, and 
NaD) to [a/Fe]. 

xiv) We compared the line-strengths of the new models 
with Galactic Globular Gluster data. We hnd good agree¬ 
ment with ]Mg/Fe] determinations obtained from individual 
stars in these clusters. The clusters with the highest metal- 
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licities fall in between the a-enhanced and the scaled-solar 
model predictions in the Mg6 - (Fe) index-index diagram. At 
very low metallicities the differences between the index val¬ 
ues obtained from the various models, i.e. base, a-enhanced 
and scaled-solar, are far less prominent. This allows us to 
conclude that the base models are also suitable for study¬ 
ing the line-strengths of such metal-poor stellar systems. 
We hnd that all the models fit the C24668 index very well, 
but not CN 2 (CNi). This is possibly related to the general 
presence of two stellar generations in GCs, the first with 
a canonical a-enhanced heavy element distribution and the 
second characterised by an anti-correlated CN-ONa abun¬ 
dance pattern. 

xv) We match reasonably well the u — g vs. g — r and r — i 
vs. g — r optical colour-colour diagrams of luminous SDSS 
ETGs simultaneously if we use appropriate combinations of 
a-enhanced models and a varying IMF. Note that this is not 
possible using base models with a standard IMF. This result 
is in good agreement with ou r previous results obtain ed with 
the MIUSGAT base models llRicciardelli et aT]|2012h . 

xvi) We obtain good fits to Early-Type Galaxy spectra 
with varying mass and a-enhancement. The SSP spectra 
are able to match the Fe- and Mg-dominated features with 
great accuracy. However, we also find that in order to re¬ 
cover the a-enhancement with this fitting technique it is 
necessary to perform the fit within a spectral range that en¬ 
compasses these a-enhancement sensitive features. Within a 
larger spectral range this information is in part washed out. 

xvii) The model grids obtained when plotting H/3o against 
Mg6, (Fe) and [MgFe] indices, which are sensitive to Mg, Fe 
and total metallicity [M/H], respectively, provide consistent 
age/metallicity solutions for the most massive galaxies. 

xviii) Using SDSS ETG spectra stacked by velocity dis¬ 
persion we find that different Balmer line indices (and also 
different definitions of the same line) provide different es¬ 
timates of the age for old stellar populations. Neither the 
a-enhanced models nor models with varying IMF are able 
to solve these discrepancies, which might be due to differ¬ 
ent effects such as a variation of abundance ratios of specific 
chemical species in the stellar atmospheres, and/or more 
complex SFHs than pure SSPs. In addition these differences 
may originate in the theoretical stellar spectra used to ob¬ 
tain the differential corrections employed to synthesise the 
a-enhanced SSP spectra. 

xix) Using ~ 5000 high quality SDSS ETG galaxy spec¬ 
tra, with S/N > 25 per A, we find a linear relation between 
a proxy for the abundance obtained using solely scaled-solar 
models, [ZMg/ZFe]ss, and the [Mg/Fe] derived when em¬ 
ploying our new models with varying abundance. We obtain 
the following relation 

[Mg/Fe] = 0.59[ZMg/ZFe]gg(^FaSTI) (®) 

Very similar relations are obtained when using the base 
model grids, either based on BaSTI or PadovaOO isochrones. 
The method consists of plotting a good age indicator (such 
as e.g., H/3 or H/Io, or any other approach that provides us a 
reliable age estimate) against Mgfe, and against (Fe), using 
a scaled-solar model grid as a reference. The difference in 
metallicity obtained from these two diagrams, i.e. the proxy 
[Zmq/Z pe], which is non-zero when the object is [Mg/Fe] en¬ 
hanced, is appropriate for estimating the ’’true” abundance, 
without making any use of a-enhanced models. 


All these models can be retrieved from the MILES web¬ 
page at http://miles.iac.es The link includes a user-friendly 
web-tool facility that, among other things, allows for the 
construction of SED models for a variety of SFHs. This in¬ 
cludes a user-defined SFH, where it is possible to combine 
bursts with varying IMFs and [a/Fe] abundance ratios. 
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